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magnetic  field,  B,  and  electrostatic  potential  are  presented  as  a function  of  altitude.  The  results 
show  the  importance  of  the  image  transport  parameter  (kLpj/ojcj)(l+i'j  2/cocj  2)-l  (where  k is  the 
wavenumber  transverse  to  B,  Lj_  is  the  transverse  dimension  of  the  cloud,  and  v-x  and  a)cj  are  the 
ionospheric  ion-neutral  collision  frequency  and  ion  cyclotron  frequency,  respectively)  in  deter- 
mining the  magnitude  of  imaging  and  aspect  angle  of  striations  with  respect  to  B (i.e.,  striations 
take  on  a parallel  component  of  wavenumber).  Perturbations  penetrate  further  down  in  the 
presence  of  the  plasma  cloud  than  in  simple  previous  studies  which  neglected  image  transport  and 
considered  the  perturbation  mapping  from  one  region  of  the  ionosphere  to  another.  Our  results 
show  that  clouds  with  smaller  conductivity  xatios  produce  image  striations  further  down  into  the 
background  E region  ionosphere  with  a more  uniform  coupling  as  a function  of  wavelength.  It 
is  also  shown  that  there  is  a slight  dependence  of  the  E region  coupling  of  the  perturbations  on 
the  level  of  solar  activity  (solar  maximum  or  minimum  conditions)  and  also  this  E region  cou- 
pling show^a slight  dependence  on  the  extent  of  F region  coupling  above  the  cloud.  Finally, 
with  our  best  estimates  for  F region  coupling,  the  growth  rates  show  negligible  short  wavelength 
damping  due  to  ionospheric  coupling  for  the  ££ /£*  = 4 case;  whereas,  for  the  ££  /Zj,  = 0.2  case 
the  ionospheric  coupling  is  significant  with  respect  to  structuring  the  spectrum  at  short  wave- 
lengths. 
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COUPLED  BARIUM  CLOUD-IONOSPHERE  SYSTEMS 
4.  STRIATION  PENETRATION  INTO  AN  INHOMOGENEOUS  IONOSPHERE 


I.  Introduction 

In  the  past  few  years  much  theoretical  and  computational  work 
has  been  devoted  towards  understanding  the  behavior  of  striations 
produced  by  barium  clouds  released  in  the  ionosphere  [Linson  and 
Workman,  1970;  Simon,  1970;  Volk  and  Haerendel,  1971;  Simon  and 
Sleeper,  1972;  Perkins  et  al»,  1973;  Zabusky  et  al.,  1973;  Lloyd 
and  Haerendel,  1973;  Shiau  and.  Simon,  1973;  Goldman  et  al.,  1974, 
Scannapieco  et  al.,  1974;  an<i  Perkins  and  Doles,  19753  * In  studying 
striation  growth  coupled  to  the  background  ionosphere  it  is  important 
to  know  how  the  induced  potential  (or  electric  field  transverse  to 
the  magnetic  field)  varies  along  the  magnetic  field  line  as  one  moves 
away  from  the  localized  region  of  the  barium  cloud.  Some  investigators 
[Far ley , 1959 > 19^0;  Spreiter  and  Briggs,  I96I;  Swift,  1972]  have 
treated  the  related  problem  of  mapping  transverse  sinusoidal,  spatially 
varying  electric  fields  produced  in  one  region  of  the  ionosphere  to 
other  regions.  These  works  characteristically  neglect  the  density 
induced  by  the  presence  of  these  electric  fields  on  the  mapping 
(imaging)  of  these  electric  fields.  Moreover,  the  mechanism  producing 
the  primary  electric  fields  is  not  considered  (a  static  problem  is 
investigated,  i.e.,  no  time  variation  of  the  fields).  Nevertheless, 
such  studies  have  shown  in  detail  [Farley,  1959*  19^0;  Spreiter  and 
Briggs,  196l3  that  for  transverse  wavelengths  ^ 1 km  the  magnetic 
field  lines  act  as  equipotentials,  i.e.,  shorter  wavelength  electric 
field  perturbations  don’t  couple  effectively  from  one  region  of  the 
ionosphere  to  another^  _ 

Note:  Manuscript  submitted  June  24,  1975. 
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Most  ideas  about  the  potential  variation  along  magnetic  field 
lines,  as  a function  of  transverse  wavelength,  are  based  on  these 
earlier  works.  Moreover,  the  multilevel,  two-dimensional  theoretical 
and  numerical  simulation  studies  of  barium  clouds  coupled  to  the  back- 
ground ionosphere  [Lloyd  and  Haerendel,  1973;  Goldman  et  al . ; 197^> 
Scannapieco  et  al.,  197^-]  have  assumed  that  the  electric  fields  at 
the  cloud  level  and  the  background  ionosphere  are  the  same.  That  is, 
the  electric  fields  have  only  transverse  spatial  dependence  in  the 
two  dimensions  perpendicular  to  the  magnetic  field  and  no  variation 
along  the  magnetic  field  lines.  Estimates  of  striation  behavior, 
based  on  linear  theory  and  including  coupling  to  the  background 
ionosphere,  have  been  done  by  Volk  and  Haerendel  [l97l]>  and  Perkins 
et  al . [ 19731 • However,  both  sets  of  authors  have  assumed  that  the 
ionosphere  is  uniform,  though  different  in  density  and  Pedersen 
conductivity  from  the  barium  cloud  (inhomogeneity  region).  Also, 
high  altitude  models  were  used,  i.e.,  the  ratio  of  ion-neutral 
collision  frequency  to  ion  cyclotron  frequency,  v/u)^  was  assumed 
small  for  both  the  cloud  and  ionosphere.  Indeed,  in  the  section  on 
"image  striations  in  a homogeneous  background"  in  Volk  and  Haerendel 
[1971],  simple  conclusions  are  drawn  based  upon  the  simplifying 
assumptions  of  homogeneous  ionosphere  and  equal  v/wq  for  the  barium 
cloud  and  the  ionosphere.  Although  these  assumptions  may  be  valid 
for  potential  modes  which  are  localized  along  the  magnetic  field 
lines  in  the  vicinity  of  the  barium  cloud,  they  are  not  in  general 
valid.  Volk  and  Haerendel  [l97l]  state  that  the  results  of  their 
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simplified  analysis  on  image  effects  on  striation  growth  show  the 
need  for  considering  a real  ionosphere.  Swift  [1972]  has  also  noted 
that  his  results  and  essentially  those  of  Farley  [ 1959 ^ 19^0]  and 
Spreiter  and  Briggs  [1961]  have  to  viewed  with  caution  when  transverse 
inhomogeneities  (e.g.,  barium  clouds)  are  present. 

With  all  this  in  mind,  the  purpose  of  the  present  study  is  to 
assess  the  interaction  between  the  earth's  ionosphere  and  striations 
which  are  generated  in  a relatively  localized  region  of  the  ionosphere 
by  density  gradients  in  a barium  cloud.  Our  study  is  posed  as  an 
eigenvalue  problem  which  while  determining  how  far  (along  the  magnetic 
field)  striations,  created  in  a barium  cloud,  will  penetrate  into  the 
background  inhomogeneous  ionosphere  also,  self-consistently,  solves 
for  their  growth  rate.  All  this  is  done  as  a function  of  transverse 
wavelength,  Pedersen  conductivity  ratio  (cloud  to  ionosphere),  and 
ambient  ionospheric  conditions  (for  real  ionospheres).  For  cases 
where  the  integrated  Pedersen  conductivity  of  the  barium  cloud, 
is  large  compared  to  the  integrated  Pedersen  conductivity  of  the 
background  ionosphere,  S^1,  the  growth  rate  of  the  instability  depends 
on  barium  cloud  parameters  and  is  essentially  known.  However,  the 
self-consistent  imaging  at  different  altitudes  in  the  ionosphere  has 
not  been  fully  treated.  For  cases  where  L ^ is  comparable  to  or 
smaller  that  £ neither  the  striation  growth  or  the  self-consistent 
imaging  has  been  fully  treated. 

We  present  in  the  second  section  the  model  and  theory;  in  the 
third  section  numerical  results  and  discussion  of  the  eigenmodes  and 
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eigenvalues  and  modifications  due  to  partial  coupling  to  the  F 
region;  and  in  the  fourth  section  summary  and  conclusions. 

II.  Model  and  Theory 

In  our  model  the  magnetic  field,  B,  is  taken  to  be  vertical 
and  in  the  z direction.  Furthermore,  the  barium  cloud  is  taken  to 
be  localized  in  z at  altitudes  high  enough  so  that  one  does  not 
have  to  be  concerned  with  the  variation  of  the  electrostatic  potential 
over  its  length  along  z.  Also,  there  is  a spatially  and  temporally 
constant  background  ambient  electric  field,  and  neutral  winds  are 
neglected. 

At  each  point  in  space  the  electron  and  ion  densities,  n, 
assumed  equal  under  the  standard  long  wavelength  (quasi-neutral) 
assumption,  evolve  respectively  according  to  the  electron  density 
equation 
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and  the  ion  density  equation. 
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where  z = B/|b|,  v denotes  collision  frequency,  denotes  cyclotron 
frequency,  e is  the  electronic  charge  (positive),  T,  is  temperature  in 
energy  units,  the  minus  subscript  denotes  electron  quantities  (so  that  v 
appearing  alone  will  mean  ion  collision  frequency  and  similarly  for 
oa^),  a is  the  Pedersen  conductivity  taken  as 
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and  \|r  is  related  to  the  ordinary  induced  electrostatic  potential 
by 
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Equations  (l)  and  (2)  are  readily  derivable  from  Eqs.  (3)  and  (4) 
of  Goldman  et  al.,  [ 197^3  and  are  in  the  E^xB  drift  frame. 

We  view  the  striations  as  due  to  the  first  order  variation  in 
a zero  order  background.  On  expanding  Eqs.  (1)  and  (2)  with 
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and  for  the  first  order  terms 
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where  we  have  used  ay  = nep/m  v . 

We  neglect  any  altitude  variation  within  the  barium  cloud  and 
assume  that  the  striations  are  localized  in  a region  perpendicular 
to  B with  barium  density  given  by  n^Q(x),  Further  we  assume 
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Eq  - = E(x).  The  assumption  bE/by  = 0 is  necessary  for  harmonic 

variation  in  the  y direction.  The  assumption  bE/bz  = 0 is  based  on 
bWbz  = 0,  which  is  reasonable  for  striations  with  variation  exp  [iky] 
such  that  k L±  » 1 where  L±  = |bln  n^/hx}  1 . That  is,  we  expect 
that  the  zero-order  potential  associated  with  L±  is  essentially 
constant  in  z over  distances  for  which  the  first  order  striation 
potential,  associated  with  a more  rapid  variation  perpendicular  to 
B,  is  varying.  The  background  ionospheric  density  is  taken  to  depend 
only  on  z and  represents  an  extension  over  the  work  of  Volk  and 
Haerendel  [l97l]  and  Perkins  et  al.  [ 19731  who  assumed  a constant 
background  ionosphere.  This  neglects  any  zeroth  order  image  clouds 
in  the  background  ionosphere  (see  Goldman  et  al.,  197^;  and  Scannapieco 
et  al . , 197^ ) and  is  certainly  true  at  early  times  in  the  barium 
cloud  development.  It  is  also  interesting  to  note  that  neglect  of 
density  gradients  perpendicular  to  B in  the  background  ionosphere 
appears  reasonable  for  the  inhomogeneous  steady  state  consistent 
with  our  assumptions  (see  appendix  A for  discussion).  The  zero-order 
background  may  in  principle  be  varying  in  time;  since  the  striation 
growth  rate  at  any  given  time  depends  only  on  the  zero-order  parameters 
at  the  same  time,  this  effect  will  be  neglected. 

With  the  preceding  discussion  in  mind  we  take 

Vt)[  =[Vi|rJ(x)  +7[  M2)elky+Yt]  (10) 
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By  omitting  the  x variation  in  first  order  quantities  we  overestimate 
V by  a fraction  proportional  to  (kLx ) 1 due  to  nonlocalization  of 
the  mode  [Goldman,  197l]  and  neglect  the  modification  of  the  insta- 
bility by  current  along  the  density  gradient  [Perkins  and  Doles, 
197b]*  The  latter  effect,  according  to  the  authors,  can  (at  least 
in  the  absence  of  diffusion)  vanish  completely  with  appropriate 
direction  of  the  electric  field  for  the  case  of  two-dimensional 
striation  variation  transverse  to  the  magnetic  field. 

Upon  substituting  Eqs.  (10)  and  (11)  in  (8)  and  (9),  going  to 

A 

a reference  frame  moving  with  v = -c  V’lr^x  .z/B,  using  kL^  > 1 and 
equivalently  < kE  (valid  for  |vi|i  |/|e|s  1),  and  expressing  nx 

in  terms  of  'll  1 (through  Eq.  (9))  as 


with  £ = E#y  + (v/io^E'X  (x  and  y are  unit  vectors  in  the  x and  y 
directions,  respectively),  we  finally  obtain 
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In  Eq.  ( 14)  we  have  that  6^  = 1 for  z^  ^ z ^ zlf  where  z^  is  the 
bottom  of  the  barium  cloud  and  zx  is  the  top  of  the  barium  cloud 


At  this  juncture  we  should  point  out  that  Eqs.  (13)  and  (l4)  outside 

the  barium  cloud,  i.e.,  6^  =0  reduce  to  the  equation  studied  by- 

Farley  [1959,  i960]  and  Spreiter  and  Briggs  [1961]  if  we  neglect 

the  factor  -ikL  (v/u)  )/(1+(v/<ju  )2)  in  the  denominator  of  Eq.  (14). 
i*  c c 

This  term  arises  from  the  contribution  of  image  density  transport 
in  Eq.  (9),  i.e.,  in  determining  the  relation  between  and  ti 
for  the  ionospheric  plasma  and  including  effects  due  to  na  (essentially 
the  terms  involving  a on  the  RHS  of  Eq . (9)). 


(i.e.,  within  the  barium  cloud)  and  6^  = 0 otherwise  (note  that 
the  subscript  b refers  to  barium)  and  also 


Pi 
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Equivalently,  if  we  neglect  the  n^  type  contributions  to  the  current 
in  the  divergence  of  the  current  equal  zero  equation,  then  we  get  the 
results  of  Farley  [1950,  i960]  and  Spreiter  and  Briggs  [ 1961] . That 
is  to  say,  we  have  from  Eq.  (13),  in  the  ionosphere, 

7-i  = r>  (16) 
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and  we  have  taken  the  perturbed  current  as 
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and  V a = iky  a , etc.  However,  for 
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i.e.,  only  contributions  from  the  perturbed  potential  i|>  t , we  have 
Eq . ( 16),  but  here 


T = k2o  1 (20) 

po 

Equations  (13)  and  ( lb ) outside  the  barium  cloud  are  equivalent  to 
equations  (6 1)  and  (62)  of  Perkins  et  al.  [ 1975]  if  we  assume  v/u)^  <<:  1 
and  is  independent  of  altitude  and  noting  that  we  have  no  x variation 
for  n2  and  ^ . Finally  Eqs.  (13)  and  (14),  in  the  background  ionosphere, 
are  equivalent  to  Eq.  (59)  °f  Volk  and  Haerendel  [l97l]  if  we  assume 
kL,Wu)cYi  <<:  1 and  we  take  a homogeneous  ionosphere. 

III.  Numerical  Results  and  Discussion 
On  defining 


K 


bz 


(21) 


we  can  rewrite  Eq.  (13)  as 


i pr  bK 
bz  bz 


(22) 


Equation  (22)  is  solved  numerically  for  the  case  of  a barium  cloud 
centered  at  an  altitude  of  200  km  and  assumed  to  extend  from  190  km 
(ZQ)  to  210  km(z1),  with  L±  = 1 km,  in  an  ionosphere  defined  for 
z ^ 100  km.  Parameters  for  the  ionosphere  are  basically  those  given 
by  Hanson  [1961]  for  nighttime  sunspot  minimum  and  maximum  (see 
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Figs.  1-3;  note:  that  the  curve  in  Hanson  [1961]  for  nighttime  sunspot 


minimum  for  the  Pedersen  conductivity  is  in  error  below  the  E region 

peak  and  in  any  case  our  Figs.  1-3  are  consistently  based  on  the 

curves  for  n^  and  v given  by  Hanson  I1961]  and  show  the  ambient 

ionosphere  that  we  used)  and  we  take  the  ionospheric  temperature,  T , 

in  our  equationsto  be  constant  so  that  by^/bz  = 0 (see  Eq.  (15))- 

^The  variation  in  y^,  due  to  the  z dependence  of  § (through  the  term 

v ( Vb  \ 

VE  /a)  ) is  only  significant  for  ( h — — - — ) ^ 1, 

x c'  j & (D  V (D  , eE  L / 9 

c v cb  y x/ 

[see  Eq.  (A-4)],  i.e.  for  z ***  100  km  where  the  Pedersen  conductivity 
is  already  small;  hence  it  will  be  neglected. | 

The  variable  K represents  the  downward  current  parallel  to  the 
magnetic  field  lines  due  to  the  striations.  Below  about  110  km 
there  is  a sharp  fall-off  in  the  Pedersen  conductivity  (see  Fig.  l); 
hence,  it  appears  reasonable  to  terminate  the  ionosphere  at  100  km. 
We  assume  no  current  flow  outside  of  the  ionosphere;  therefore 
K(°°)  = K(100  km)  = 0.  The  condition  K(°°)  = 0 coupled  with  estimates 
of  potential  variation  for  the  ionosphere  based  on  decreasing  values 
for  the  electron-neutral  and  ion-neutral  collision  frequencies  at 
high  altitudes  suggest  that  we  treat  the  region  above  the  barium 
cloud  as  being  at  constant  (unit)  potential,  i.e.,  ^i(z)  = 1 for 
z ^ 190  km(zQ).  This  assumption  is  investigated  numerically  as 
concerns  its  effects  on  growth  rates  and  modal  structures  (see  cases 
IV  and  V in  this  section). 

The  boundary  conditions  on  Eq.  (13)  may  be  written  as 


* 
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00 


(23a) 


K(zq)  = -J  r(z'  )dz 


Z 

O 


K(  100  km)  = 0 


(23b) 


(23c) 


(remembering  that  zq  = 190  km)  where  (23a)  and  (23b)  are  the  current 
conditions  at  infinity  and  100  km,  respectively  and  (23c)  shows  that 
the  potential  is  taken  to  be  unity  at  190  km.  Since  Eq.  (22)  is 
second  order  and  there  are  three  boundary  conditions,  it  is  clear 
that  Eqs.  (22)  and  (23)  pose  an  eigenvalue  problem.  This  problem 
has  been  solved  using  the  B0DEL  program  (see  Appendix  B)  with  129 
points  in  the  range  100  km  ^ z 5 190  km. 

Results  for  perturbed  potential,  number  density,  and  parallel 
velocity  (all  as  a function  of  altitude)  and  growth  rates  for  wave- 
lengths from  0.2  km  to  1.0  km  (in  steps  of  0.1  km)  and  for  wavelengths 
from  1 km  to  6 km  (in  steps  of  1 km)  are  available  for  the  models 
listed  below.  Figures  4-1 6 display  these  results  for  selected  wave- 
lengths and  different  models.  The  models  used,  hereafter  to  be 
referred  to  by  their  Roman  numerals,  are  : (i)  the  canonical  situation 
of  nighttime  solarmin  (sunspot  minimum)  conditions  [according  to 
Hanson,  1961],  a cloud  with  ratio  of  height  integrated  Pedersen 
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conductivity  to  that  of  the  background  ionosphere,  E k/JJ  ^ equal 

to  4 and  full  coupling  to  the  F region  (^i(z)  = 1,  z > 190  km); 

( II ) the  canonical  situation  but  with  E ^/E  * = 0.2;  (III)  the  canonical 

P p 

situation  but  with  the  value  zero  substituted  for  - ikL ,v/u)  /(l+(v/u>  )2) 

± c c 

in  the  denominator  of  Eq.  (14),  so  that  this  effectively  suppresses 

the  contribution  of  image  density  transport  in  determining  the 

relationship  between  nx  and  for  the  ionospheric  plasma  and  provides 

the  limit,  for  the  ionosphere,  appropriate  to  the  earlier  calculations 

of  potential  variation  along  field  lines  [Farley,  1959,  19&0;  Spreiter 

and  Briggs,  I96I;  and  Swift,  1972];  (IV)  the  canonical  situation  but 

with  solarmax  conditions  L Hanson,  1961] ; and  (V)  the  same  situation 

as  (IV)  but  with  the  assumption  ti(z)  = 1 for  190  km  < z < 210  km 

and  ti(z)  = 0 for  z > 210  km.  The  difference  between  (IV)  and  (V ) 

is  presented  to  provide  an  upper  bound  on  the  uncertainty  due  to  F 

region  coupling.  It  is  worthwhile  to  point  out  at  this  juncture 

that  for  solarmin,  from  Fig.  1,  that  E^1  =0.19  mhos  while 

E^  (above  the  barium  cloud)  = 0.04  mhos;  whereas,  for  solarmax 

Ep  = 0.75  mhos  while  E * (above  cloud)  =0.58  mhos. 

Eigenmodes  • ^ is  normalized  so  that  ty1(190)  = 1;  whereas  n^  and  K are 

presented  in  units  of  (k'lti/^).  For  E^/E^  » 1 one  expects 

Wx/g  ~ 1,  and  for  Z?3 /E*  « 1 one  expects  ki^/^  ~ (S^/E1 ) , before 

P P P P 

non-linear  effects  on  modal  structure  set  in. 

* The  phase  variation  present  in  all  cases  except  III  (see 
Figs.  10-14)  has  the  following  interpretation:  A modal  structure 
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for  ih  (or  v 

iz 


— — or  Hi  ) can  be  written  as 
N e x ' 
o 


(tir  + i ♦li)  exp  (iky  + Yt) 
or  equivalently  as 

| ilij  | exp[  iky  + i§(z)  + yt] 


with  | til  = Vtlr^  + and  sin  §(z)  = txi/|  til  , cos§(z)  = tlr/|  tj  . 

Hence  the  phase  ( = [ky  + §(z)  + is  the  same  at  a given  time  for 

locations  y(190  ) and  y(z),  just  below  I90  km  and  at  z,  respectively, 
provided: 


y(z)  = y(l90")  + 


CJ.(1SP'.)  - X + 

2tt 


where  m is  integral  or  zero  and  X = 2n/k  is  the  wavelength  perpendicular 
to  the  magnetic  field  in  the  y direction;  i.e.  the  harmonic  modal 
structures  shift  along  they-axis  as  z varies.  This  variation  is 
shown  for  case  I at  .8  km  in  Fig.  (ba)  and  for  case  II  at  .4  km  in 
Fig.  (4b).  (The  shifts  are  reflected  around  y = o for  a negative  B 
field  as  is  appropriate  for  the  northern  hemisphere.)  It  is  interesting 
to  point  out  that  in  both  cases  presented  in  Fig.  4 the  density 
variation,  nly  just  under  the  cloud  (190  km)  is  out  of  phase  with  the 
striations  in  the  cloud  by  approximately  90°  (1/4  of  a wavelength). 
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This  is  in  agreement  with  the  homogeneous  ionosphere  theory  given  by 
Volk  and  Haerendel  [l97l]*  However,  the  subsequent  density  variation, 
as  one  descends  in  altitude,  is  not  as  pronounced  nor  amenable  to 
as  simple  an  interpretation  as  given  by  their  Fig.  8 (or  the  text 
pertaining  to  that  figure),  i.e.,  real  ionospheres  are  somewhat  more 
complicated  than  homogeneous  ionospheres  in  the  treatment  of  the 
image  problem. 

One  notes  that  one  may  define  a wave-vector 

k ( z ) = k ( z ) + ik  . ( z ) 
z 7 zr  zi 

such  that 


kzr(z ) = d§(z)/dz 


and 


kzi.(Z)  = ‘ d lnl  ( z ) I /dz 

to  characterize  the  z-variation  of  ih  (or  v or  nO.  With  a 

1 z 

spatially  homogeneous  ionosphere,  dkz/dz  = 0 (see  for  example  Volk 

and  Haerendel,  1971)  and  so  one  has  k and  k . being  constant  as  a 

function  of  altitude  (although  not  equal).  In  our  study  k and 

zr 

k^  are  functions  of  altitude  and  in  general  are  not  simple  (this 
can  be  seen  by  examining  Figs.  5- 1*0* 
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For  Case  I the  potential  is  progressively  less  localized  with 
increasing  wavelength  (see  Fig.  5a)-  At  all  wavelengths  there  is 
an  imaginary  component  to  the  potential  (see  Fig.  10a);  at  relatively 
shorter  wavelengths  this  component,  which  can  dominate  over  the  real 
component  of  the  potential  at  the  same  altitude,  can  be  as  large  as 
0.3  of  the  value  of  the  potential  at  I90  km.  The  values  for 
vlz( 190 ) which  are  typically  3*5  x 10b  cm/sec  for  k = 1,  are 
closely  spaced  for  \ ^ 1.0  km  and  decrease  in  magnitude  as  X increases 
above  1.0  km  (see  Fig.  5*0  The  imaginary  components  of  vlz(190)  are 
roughly  2/3  of  the  real  components  of  vlz(l90)  (see  Fig.  10b).  It 
can  be  shown  that 


(♦1  = 0 = 


n ev 
o 


lz 


(k*i(i90)/£  = l) 


(24) 


where  K(^i  = l)  is  the  value  of  K appropriate  to  1^(190)  = 1 and 

viz(k^1(190)/<f  = 1)  is  the  value  of  vlz  appropriate  to  k^  ( I90  )/£  = 1, 

i.e.  the  value  0 f vlz  plotted  in  our  figures.  Since  the  contribution 

of  the  region  below  the  barium  cloud  to  Yi  occurs  solely  through 

K(ti  = l),  see  Eq.  (28),  one  sees  that  one  can  parametrize  the  entire 

region  below  the  barium  cloud  by  an  effective  complex  Pedersen 
^ /_  Kr  Ki  \ 

conductivity,  E l = — + i 1,  which  varies  roughly  as  k 1 for 

X ^ 1.0  km.  By  Eq.  (24)  the  plots  of  vlz,  particularly  vlz(190). 


directly  yield  the  values  of  E, and  the  locations  of  the  digits,  m, 
m < ()  give  effective  integrated  conductivities  between  100  km  and 


1Y 


the  location  denoted  by  (100  + 10  m)  km.  The  values  for  n* , expressed 
analytically  in  terms  of  by  Eq.  (12),  show  that  the  imaging  for 
relatively  long  wavelength  disturbances  peaks  at  the  Pedersen  con- 
ductivity maximum,  between  z = 120  and  z = 130  km,  but  that  the  imaging 
from  shorter  wavelengths  peaks  at  higher  altitudes  and  is  slightly 
weaker  (see  Figs.  5C  and  10c). 

The  results  from  case  II  for  the  potential,  i|t]  , are  strikingly 
less  localized  than  those  from  case  I (see  Figs.  6a  and  11a).  This 
is  consistent  with  the  lower  value  of  | Vi  I (see  Fig.  15)  which  leads 
to  decreased  values  of  and  ^ence  l through  Eq.  (13)]  to 

larger  scale  lengths  for  the  variation  of  i|i  in  the  z-direction,  as 
well  as  the  relatively  larger  size  of  V-^ 1 Ya  | for  Case  II  [see 
Eqs.  (24a)  and  (24b)  and  subsequent  discussion].  The  behavior  at 
.2  km  for  Case  II  is  particularly  interesting  in  view  of  the  non- 
monotonicity of  i|ii . Such  a result  is  possible  with  Eq.  (13)  if 
there  exist  values  of  z where  Pvii  > | Y1 | (with 
P = kLxW[u)c ( 1 + v2/(uc2)]).  To  see  this  we  note  on  multiplying 
Eq.  (13)  by  'll*,  integrating  between  100  km  and  z,  using  a o6'l>i/bz(l00)  = 0, 
performing  the  analogous  operations  on  the  complex  conjugate  of  Eq . 

(13)  and  adding  the  results  that  one  obtains: 


(24a) 
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where  for  z < z : 
o 


(r  + r*)/k2c 


po 


2|yJ2  - 2YUP 
I Yih  + - SY^P 


(24b) 


In  the  zero  temperature  limit  when  the  left  hand  side  of  (24 a) 

/*  Z 2 

is  equal  to  Re  2 I ji#Eidz'  the  condition  a||0^l,lril  /&z<  0 corresponds 

^z  = 100 

to  the  existence  of  a physical  volume  100  ^ z ^ z within  which  the 
ionospheric  particles  do  work  on  the  wave  rather  than  the  opposite. 

On  the  other  hand  if  one  neglects  the  contribution  in  (3  to  F from 
(24a)  and  (24b)  one  has  the  more  usual  result: 


°Ho  t>Ua|a/0*>°  • 


Further  one  notes  through  (24a)  and  (24b)  the  possibility  of  connection 

between  > Iyi|2  and  decreased  values  for  |b^/6z|  and  b|i|fi|2/6z, 

and  hence  the  possibility  of  less  localization  of  the  potential  with 

Yii  > 0,  even  if  b|^i|2/bz  > 0 for  all  z(<  zq). 

The  behavior  for  vlz  (and  hence  £)  is  markedly  different  from 

that  of  Case  I.  |vlz(190)|  is  monotonically  decreasing  as  a function 

of  increasing  wavelength  (see  Fig.  6b);  further  the  real  part  of 

v (190)  is  relatively  constant,  whereas  the  imaginary  part  of 
1 z 

vi z(190)  decreases  with  increasing  wavelength  (see  Fig.  lib).  This 

corresponds  to  a picture  of  £ whose  real  part  varies  roughly  as  k 1 

but  whose  imaginary  part  is  roughly  constant  at  short  wavelengths  and 

at  any  rate  varies  less  rapidly  with  k.  In  keeping  with  the  non- 

b i 


linear  results  for  barium  clouds  with  £ /£ 

P P 


= e « 1 [Goldman  et  al.. 
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197*0  we  would  expect  the  maximum  attainable  velocities  in  the  linear 

regime  to  be  multiplied  by  a factor  less  than  unity.  (For  e « 1, 

l/2 

this  factor  would  be  ~ e . ) The  densities  are  less  sharply  peaked 
than  those  for  Case  I and  larger  in  terms  of  (k^i/^*)  especially  at 
short  wavelengths. 

The  modal  structure  curves  for  tyj , vlz  and  nx  of  Case  III  (see 

Fig.  7)  are  universal  in  the  sense  that  they  are  independent  of  Yi 

(or  equivalent  of  the  ratio  £ ^/£  ^.)  Hence  Case  III  can  be  directly 

P P 

compared  with  both  Cases  I and  II  as  concerns  modal  structure. 

Agreement  is  closer  between  Cases  III  and  I than  between  III  and  II. 

The  results  from  Case  III  for  ifx  (see  Fig.  7a)  are  more  localized 
than  those  from  Case  I although  the  values  of  Yi  are  very  close. 
Comparison  of  Eqs.  (17)  and  (20)  results  in  the  modification  of  the 
first  term  of  Eq.  (13)  by  a factor  of  Yi/(Yi  - i-P ) in  passing  from 
Case  III  to  Case  I.  Hence  at  wavelengths  with  kL^  » 1 with  locations 
in  z with  {3  » |yi|,  b/bz (aj  <btyx/bz ) is  clearly  smaller  for  Case  I 
than  for  Case  III  and  the  potential  does  not  vary  as  rapidly  in  z. 

The  values  of  |vlz(190)|  for  Case  III  (see  Fig.  7b)  are  charac- 
teristically 20-30$  higher  than  those  for  Case  I indicating  corresponding 
larger  values  for  |L|  and  hence  lower  growth  rates.  The  constancy  of 
|viz(190)|  for  X s 1.0  km  indicates  that  |L|  varies  as  k 1 in  this 
range  of  X.  The  contrast  with  L from  Case  II  is  striking.  The 
values  of  | nx | for  Case  III  (see  Fig.  7)  are  higher  in  maximum  value 
and  more  localized  in  space  than  those  for  Case  I;  ni(z)  bears  a 
fixed  phase  relation  to  i|ii(z),  independent  of  z (see  Fig.  12). 
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The  results  from  Case  IV  for  ^ (see  Figs,  8a  and  13a)  indicate 
more  localization  in  the  vicinity  of  the  barium  cloud  than  for  Case  I. 
Mathematically,  for  two  ionospheres,  a and  b,  with  neutral  densities 
n such  that  n (zj/n^z)  = ot  the  modal  structure  is  the  same  for  wave- 
lengths X^  and  such  that  X^  = Q'X^,  provided  each  of  the  modes  is 
localized  in  a region  with  v/u)^  « 1.  Since  the  neutral  densities  in 
the  vicinity  of  the  barium  cloud  are  roughly  a factor  of  two  larger 
for  solarmaximum  than  for  solarminimum  conditions,  the  rough  equivalence 
for  potential  variation  in  the  vicinity  of  the  barium  cloud  for  wave- 
lengths such  that  X(Case  IV)/X(Case  1)^2  is  plausible.  At  lower 
altitudes  where  the  ionospheres  of  Cases  I and  IV  are  more  similar 
the  fall  off  in  ifx  with  decreasing  z at  the  same  wavelength  is  closer. 
This  gives  the  curves  of  |^x|  for  Case  I a more  bowed  appearance  than 
those  of  Case  IV. 

The  values  of  |v12(190)|  for  both  cases  are  similar  (see  Figs.  5b 
and  8b).  At  relatively  small  values  of  X,  v1Z£( 190 )/vlzr ( 190  ) is 
larger  for  Case  IV  (see  Fig.  13b)  than  for  Case  I (see  Fig.  10b). 

This  effect  which  also  is  pronounced  for  short  wavelengths  in  Case  II 
is  apparently  related  to  relatively  large  values  of  y3^,  as  can  be 
seen  from  the  analytic  expression: 

Yir2  + Yxi(Vii-P)  + iPYir 

Yir2  + (Y^-P)2 

(25) 

♦i(Oa  (Od*'  , 


ViZ'190^  n (190)< 


*100 
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which  indicates  the  possibility  of  negative  contributions  to  vlzr 
for  regions  of  space  such  that  py  > I Y, I • At  short  wavelengths 

ii  •*- 

the  maxima  for  | nx | in  Case  IV  (see  Fig.  8c)  are  larger  and  occur 
at  higher  altitudes  than  in  Case  I (see  Fig.  5°)*  This  corresponds 
to  the  larger  Pedersen  conductivity  under  solarmaximum  conditions 
at  higher  altitudes  in  the  range  100  - z ^ lyO  as  well  as  the  greater 
localization  of  i|i  . At  large  wavelengths  the  agreement  is  closer 
with  the  maxima  in  | n^ | still  at  slightly  higher  altitudes. 

The  results  from  Case  V for  ^ (see  Figs.  cJa  and  l4a)  are  very 
similar  to  those  from  Case  IV  except  at  short  wavelengths  where  there 
is  more  localization  of  the  potential  near  the  barium  cloud.  At  short 
wavelengths  one  notes  that  y^  from  Case  IV  is  relatively  much  larger 
than  y ^ for  Case  V (see  Fig.  15)*  By  Eqs.  (24a)  and  (24b)  this 
suggests  a less  rapid  decrease  of  Itxl  with  increasing  distance  from 
the  barium  cloud.  The  larger  values  of  ^v  (190)|  at  short  wavelengths 
in  Case  IV  (see  Figs.  8b  and  9b)  are  attributable  to  the  larger 
magnitudes  of  vlz^(190)  which  appear  to  be  due  to  the  large  values  of 
for  z ^ I50  km.  The  larger  values  of  n3  (see  Figs.  8c  and  9C)  at 
short  wavelengths  in  Case  IV  are  attributable  to  the  less  rapid  fall- 
off  in  i|r  away  from  z = 190  as  well  as  the  more  negative  values  of 
the  imaginary  part  of  y. 

Coupling  above  the  barium  cloud  and  eigenvalues.  In  this  section,  we 
first  present  conclusions  obtained  in  Appendix  D regarding  electrostatic 
coupling  above  the  barium  cloud.  We  then  apply  these  notions  to  the 
estimation  of  induced  velocities  parallel  to  the  magnetic  field, 
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density  inhomogeneities , and  the  calculation  of  eigenvalues  taking 
partial  coupling  to  the  F region  conductivity  maximum  into  account. 
The  reader  is  referred  to  Appendix  D for  analytic  details. 

For  z >zly  the  treatment  of  Appendix  D indicates  for  X ^ 0.2  km 
and  the  barium  clouds  of  Cases  I and  II,  or  for  X ^ 0.5  km  and  the 
barium  clouds  of  Cases  IV  and  V,  that  the  dominant  part  of  ^i(z)  is 
given  by 


iliio(z)  •*  Mzl) 

1 + or  ( zx  ) 


(26) 


This  result  is  in  close  agreement  with  the  physical  model  in 
which  \Jf joo , the  potential  [assumed  constant  and  equal  to  at 

the  F region  conductivity  peak  (z  > z*),  causes  an  integrated  horizontal 

c 60 

current  divergence  -/  r\|rl00dz  for  z*  ^ z ^ 00  which  in  turn  draws  a 

Jz* 

r 

vertical  current  -/  dz  into  the  region  z > z*  along  the  magnetic 

l z* 

field  lines.  If  there  is  no  significant  perpendicular  shorting  of 
current  between  zx  and  z**,  so  that: 


00 


r^oodz 


dz* 


one  has: 


-o 


llo 


d 

di~ 


j 
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from  which  (26)  follows,  provided  one  takes  kL^v/cuc  « 1 for  z > z3  . 

For  Cases  I and  II,  we  have  ot(zx)  = 3.2  x 10  2,  and  for 

Cases  IV  and  V,  ct(zx)  = X 2,  where  X is  in  km.  Further  we  note 
a(z)  <a(Zi)  for  z > zX}  so  that  d\|f10/dz  < 0,  z > zx . For  the  F 
region  conductivity  peak  region  [z  > z*  (=  280  km  at  solarmax), 
z > z* ( = 270  km  at  solarmin),  see  Fig.  l3,  we  have  or(z)  ^ 0.2.  Hence 
for  solarmin  conditions  a transition  from  little  coupling  to  the  F 
region  conductivity  maximum  to  full  coupling  is  indicated  as  X increases 
through  0.2  km;  while  for  solarmax  conditions  the  transition  is  indicated 
for  X « 1.0  km. 

In  analogy  to  Eq.  (25)  one  has  on  neglect  of  contributions  in  P 
which  are  expected  to  be  small  by  virtue  of  the  discussion  of  Appendix 
D: 

2 C 00 

Viz(zi)  = - en0(zx)  J Vl(Z')dZ'  ' (27a) 


which  by  (28'.  and  (2fiP  ),  q.v.,  yields  [with  k<lii(z1)  = £]: 


Vlz(Zl) 


k SI  a dz  ' 

yzi  po 

eno(Zi)  [l  + a(zx  )] 


(27b) 


Within  the  range  of  validity  of  Appendix  D under  solarmin  conditions, 
the  maximum  value  for  Viz(z!),  z^  = 210,  occurs  for  X = 0.2  km  where 
vlz(zi)  = 5 x 105  cm/sec;  under  solarmax  conditions  the  maximum  value 
occurs  for  X * 1.0  km  where  viz(zx)  = 1.4  x 10f  cm/sec. 
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From  Eq.  (12)  on  neglecting  terms  of  order  kL  v/co  in  the 

J.  c 

denominator,  one  has  that 

ni(z)  « - k2apQi)i1(z)/ev.  (27c) 

For  all  cases  covered  by  Eq.  (26),  n^z)  is  largest  in  magnitude  for 

f\j  rsj 

Z > z*.  On  using  y = y c£/BLx  with  y of  order  unity,  Eq.  (27c)  leads 
to  the  result 


ni(z)  - [ 1 + g(z )]  k^Zi}  kL  _v 
nQ  Ll  + «(zi )]  X cuc 


(27d) 


Hence  for  striations  with  kijti  (zx  )l£y  **  1,  as  is  expected  to  be  the 
case  at  the  limit  of  a linear  treatment,  nj.(z)/n  can  be  as  large  as 
3 x 10  2,  (at  X <=*  1 km  and  z «=*  300  km.  under  solarmax  conditions, or 
at  X 0 .2  km  and  z « 280  km  under  solarmin  conditions . ) 

First,  as  concerns  eigenvalues,  it  is  interesting  to  note  that 
the  eigenvalue  agreement  between  Cases  III  and  I (see  Fig.  15)  is 
considerably  closer  than  between  Case  III*  (Ij^ = 0.2,  but  neglect 
of  image  density  transport)  and  Case  II  (see  Fig.  15a).  In  comparing 
these  cases  (in  Fig.  15a  and  remembering  that  III  and  m*  neglect 
the  effect  of  images  back  on  striation  growth,  i.e.,  no  image  trans- 
port) we  do  not  find,  for  a real  ionosphere,  the  pronounced  effect  of 
image  striations  reducing  the  growth  rate  that  Volk  and  Haerendel 
[l97l]  found  for  the  homogeneous  ionosphere  (see  their  Fig.  9)*  f°r 
Case  IV  (full  coupling  to  the  F region),  as  depicted  in  Fig.  15a,  for 
long  wavelengths  - vlr  =0.8  which  represents  the  short-circuited 
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growth  rate  since  Ep  /E^*"  = 4.0.  For  Case  V,  (no  coupling  to  the  F 

region)  Fig.  15a  shows  that  for  long  wavelengths  -y  0.95*  This 

b i. 

can  be  understood  by  the  fact  that  although  we  chose  Ep  /Ep  =4.0 
in  this  case,  there  is  no  coupling  to  the  F region  where  most  (" 
of  the  integrated  conductivity  occurs  (recall  that  for  solarmax 
2^*  = 0.73  mhos,  and  Ep*  (above  cloud)  = O.58  mhos).  Consequently, 
the  effective  integrated  conductivity  of  the  ionosphere  for  long  wave- 
lengths is  0.15  mhos  so  that  effectively  2 i ^ w 0.05  and  this 
results  in  a short-circuited  value  of  0.95  for  -Yir* 

To  obtain  modifications  on  yx  and  y due  to  partial  coupling  to 
the  F region,  we  note  that  Eq.  (13)  yields,  on  integrating  upwards  to 
"con  from  z the  bottom  of  the  barium  cloud,  using  Eq.  (21)  and 

neglecting  6,  v.  /uj  , kL 
b b cb  1 


sp  (yx  + 1) 

yx  - ikLxvb/a)c~  + 

1 + V 2/(l)2  , 
b cb 
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0Yiti(z')dz' 


K. 


Yx 


ikL.  v/o) 
x c 

1 + V2/(l)  2 

c 


4-  + 


iK. 

1 


= 0,  , (28) 


with  K(k,ZQ)  = K^+ iK^  and  t|rx(zo)  = 1.  From  Eq.  (26)  and  the  dis- 
cussion after  Eq.  (D-l4b)  one  sees  that  the  bulk  of  the  contribution 
to  the  integral  of  Eq.  (28)  comes  from  the  region  with  z > z*  where 
kL±  v/Uc(l  + v2/(i)c2)]  « yx  and  t|»x(z')  « 'lr1o(“)  ^ 1/Cl  + <*(zx)]. 
Hence  the  integral  may  be  approximated  by 


2b 


(28') 


/> 


tlo(°°)  dz' 


which  in  turn  may  be  approximated  by  Lp  /[  1 + a ( zx )]  where 


f _ C i . * 
= / o dz  . 

> Jz  p° 


Then  Eq.  (28)  yields 


Vl"l 


-<  E b[E  b + E i(k)]  + 

' P P Pr 


kL. 


^ E [ (k)E  b 

d + V/(°cb)  u,cbP  P’ 


(28  ' ) 


(29) 


+i  (e  ^ *(k)  + |e  *(k)[L  b + s *(k)] 

V p pi  (1  + vb2/0)J)  ®cb(  pr  P pr 


where  A = [Lpb  + EprX(k)]2  + [Ep*(k)]2, 


Ep^(k)  - Kr/k2  + Spf/Cl  + cr(zi)],  and 


Ep^(k)  = K^/k2.  Eq.  (29)  and  the  resulting  values  of  y have  been 

evaluated  using  values  for  K and  K . obtained  from  Cases  1,  II  and 

Kr  Ki 

We  note  from  Eq.  (2k)  that  the  values  of  and  —^z  are  directly 


27 


related  to  viZ,  evaluated  at  z = 190  km;  further  that  vlz  evaluated 
at  190  km  is  given  in  Figs.  10b  - l4b. 

For  either  the  solarmax  or  the  solarmin  case  with  partial  coupling 
and  E 1 = 4.0  we  note  that  both  the  values  of  Yi  and  the  ionospheres 
below  the  barium  clouds  are  similar;  further  from  comparison  of  Cases  IV 
and  V at  the  same  k,  one  has  |ak|/|k|  w |Ayx|/|yi| • (Here  Ayx  specifies 
the  variation  of  an  eigenvalue  and  AK  is  the  corresponding  variation 
of  K. ) In  addition  we  have  AK  w K*(-1)Ayi,  where  K = K(Yi ) and 
K*(-l)  = bK/&Yi  evaluated  at  Yi  = -1,  so  that  with  |yi|  ~ 1, 

K'(-l)  ~ K(-l). 

On  taking  large  compared  to  other  terms  in  Eq.  (29)  as  is 

generally  appropriate  for  = 4.0,  we  may  write  Eq.  (29)  to 

b 

lowest  order  in  (L  ) as: 

P 


Yx  + 1 = 


K L kL, 

r . P _ -i 

kJ’  T+aJzT)  ~ ( 1+v^/oj^ ) 


+ i 
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cb 


Ki 

“cb  W 
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— r»  + , r 

Tc^  1-tOf  ( ZX  ) 


(29a) 


with  K = K(Yi)  ^ K(-l)  + (Yi  + 1)  K*(-l).  Since  the  quantity 

(Yi  + l)  is  small  if  is  relatively  large,  the  terms  linear  in 
K#(-l),  ~ (yi  + l)2,  are  of  second  order  in  the  smallness  parameter, 
(Yi  + 1).  Hence  the  values  of  Kr  and  obtained  from  full  coupling 

can  be  used  with  partial  coupling  to  first  order  accuracy  in  (yx  + l). 
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For  the  solarmax  situation,  the  partial  coupling  calculations  for 
X ^ 0.5  km  (see  Fig.  1 6)  attach  smoothly  to  the  uncoupled  calculations 
of  Case  V (see  Fig.  15)  for  X < 0.5  km.  For  the  solarmin  situation 
the  partial  coupling  results  for  X ^ 0.2  km  (see  Fig.  16)  go  smoothly 
into  the  full  coupling  results  at  about  X = 0.8  km  (see  Fig.  15).  For 
both  solarmin  and  solarmax  situations  one  concludes  that  there  is  no 
cut-off  in  the  striation  growth  rates  at  short  wavelengths  due  to 
interactions  with  the  ambient  inhomogeneous  ionospheric  plasma. 

For  Case  II,  the  eigenvalue  estimates  are  not  as  reliable  since 
we  do  not  have  the  data  comparable  to  that  of  Cases  IV  and  V,  and  since 
the  ordering  in  (yi  + l)  does  not  obtain.  For  X £ 0.4  km  the  eigen- 
values with  full  and  partial  coupling  agree  to  within  20$,  (see  Figs. 

15  and  l6),  so  that  there  is  presumably  little  uncertainty  in  this 
range.  One  can  obtain  an  upper  bound  on  growth  rates  at  shorter  wave- 
lengths under  the  assumption  that  coupling  below  the  barium  cloud 
reduces  Y • Then  choosing  K = 0 for  X = 0.2  km  with  partial  F region 
coupling  yields  Yr  < 0*23,  which  appears  sufficient  to  establish  the 
dip  below  Yr  ^ 0.30  at  short  wavelengths. 

In  addition,  we  find,  at  wavelengths  sufficiently  short  that  the 
perturbations  can  be  regarded  as  just  localized  within  a (substantially) 
constant  region  of  the  ionosphere,  that  coupling  to  the  ionosphere 
yields  the  growth  rate: 


po 
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with  k = JjT  mc^/L±  and  apQ  'the  approximate  local  Pedersen  conductivity 

at  the  level  of  the  barium  cloud,  provided  CTpQ/apg  <<:  (see  Appendix 

E for  details).  For  solarmax  conditions,  with  £ * = 4.0,  we  obtain 

Yr  = O.98  ciJ’/BL^  at  X = 1.1  x lO^;  for  solarmin  conditions,  with 

S ^/£  * = 4.0  we  obtain  y = 0.91  c<?/BL  for  X = 47m,  and  with 
p p r x 

£p^/£p^  = 0.2  we  obtain  Yr  = - 0.9  C(f?/BL±at  X = 47m,  where  in  obtaining 
the  above  values  we  have  used  £=  5 mV/m  and  T = 103  °K. 

IV.  Summary  and  Conclusions 

In  this  work  we  have  attempted  to  provide  a critical  quantitative 

basis  for  the  study  of  the  structuring  and  growth  rates  of  striations 

associated  with  barium  clouds  in  the  ionosphere  (i.e.,  a real  ionosphere 

whose  properties  vary  with  altitude).  We  have  been  concerned  with 

three  zero-order  physical  situations  for  a barium  cloud  between  190 

and  210  km  (£  ^/£  ^ = 4.0,  solarmin  conditions;  £ ^/£  * = 0.2, 

P P ' P P 

solarmin  conditions;  £ ^/£  * = 4.0,  solarmax  conditions)  which  have 

P P 

resulted  in  five  models  (see  Section  III,  Cases  I-V ) for  striation 
behavior  below  the  barium  cloud.  These  five  models  with  specified 
behavior  above  the  barium  cloud  have  been  investigated  numerically  as 
concerns  modal  structure  and  growth  rates  for  the  region  below  the 
barium  cloud  (Section  III).  Analytic  estimates  have  then  been  used 
to  obtain  a more  precise  specification  of  the  real  (rather  than  model) 
behavior  above  the  barium  cloud  (Section  III,  part  under  eigenvalues). 

These  in  turn  have  been  used  to  provide  corrections  to  the  growth 

— _ »— — — ■*'  ' 

rates  (Section  III,  second  part)  appropriate  to  the  real  physical 
situations  approximated  by  the  original  models. 
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Comparison  of  Cases  I (Z^/E^*  = 4.0,  solarmin)  and  II 

(E  * - 0.2, solarmin)  establishes  the  dependence  of  image  structure 

on  Spk/S  1 for  a given  wavelength.  The  smaller  conductivity  cloud 

produces  image  striations  farther  down  into  the  background  E region 

ionosphere  (this  is  in  qualitative  agreement  with  the  results  of  image 

striations  in  a homogeneous  background  given  by  Fig.  9 in  Volk  and 

Haerendel,  1971)  with  a more  uniform  effective  conductivity  as  a function 

of  wavelength.  Comparison  of  Cases  I and  II  with  Case  III  (the  limit 

of  earlier  calculations  suppressing  image  density  transport;  see 

Farley,  1959,  19^0 5 and  Spreiter  and  Briggs,  1961 ) shows  that  striations 

penetrate  further  down  into  the  ionosphere  than  expected  on  the  basis 

of  Case  III  and  that  they  make  an  aspect  angle  with  respect  to  B which 

varies  as  a function  of  altitude  (variable  parallel  component  of  wave- 

h i 

number.)  Comparison  of  Cases  I and  IV  (E^  /Z^  = 4.0,  solarmax)  establishes 

the  relatively  slight  dependence  of  E region  coupling  on  the  level  of 
solar  disturbance.  Comparison  of  Cases  IV  (totally  coupled  F region) 
and  V (no  coupling  to  F region)  established  the  relatively  slight 
dependence  of  E region  coupling  on  F region  coupling  even  under  solar- 
maximum  conditions  when  the  F region  conductivity  is  greatest. 

Some  comments  with  respect  to  the  extension  of  the  modal  structure 


results  to  other  situations  are  in  order: 


1.  One  expects  that  penetration  below  the  barium  cloud  should 
increase  with  increasing  L^.  Put  another  way,  imaging  with  Case  III 
can  be  viewed  as  behavior  in  the  limit  ■*  0.  Hence  reducing  Lx 
brings  results  under  solarminimum  conditions  closer  to  those  of  Case  III. 
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2.  Imaging  with  very  large  clouds  presumably  obtains  with 

Ylr  ^ -1,  ^ 0.  Case  V offers  the  best  approximation  to  this 

situation.  The  transition  from  Case  IV  to  Case  V indicates  how  the 
limit  » i ig  approached. 

3.  As  an  extension  of  the  comparison  of  Cases  I and  IV,  for 
very  large  clouds  so  that  Ylr  ^ -1,  Y^  w 0,  and  short  wavelengths 
such  that  ^ ^ 0 for  v/a»c  « 1,  changing  the  cloud  altitude  so  that 
the  neutral  density  in  the  vicinity  of  the  cloud  is  multiplied  by  a 
factor  of  a should  result  in  the  same  behavior  below  the  barium  cloud 
for  a new  wavelength  equal  to  the  original  wavelength  multiplied  by 

a factor  of  a (provided  the  scale  height  variation  as  a function  of 
altitude  is  neglected.) 

4.  If  all  ionized  densities  are  multiplied  by  a constant  factor 
the  results  are  unchanged  (i)^,  ni/nQ,  vJZ  and  y are  the  same)  to  the 
extent  that  Coulomb  collisions  can  be  neglected  compared  to  charged 
particle-neutral  collisions. 

Analytic  treatment  for  z ^ z^  indicates  for  X ^ 0.2  km  and  solar - 
min  conditions  and  for  X ^ 0.5  km  and  solarmax  conditions  that: 

a.  is  predominantly  real  and  decreasing  with  increasing  z, 

b.  the  maximum  of  vlz  is  close  to,  but  above  zx, 

c.  the  maximum  of  nx  is  close  to,  but  above  z = z*  (at  about 
300  km ) . 

The  estimate  of  a value  of  v Z(z  ) as  high  as  1.4  x 10s  cm/sec  is 
suggestive  of  the  possibility  that  striations  may  drive  shorter  wave- 
length ion-cyclotron  modes  [Kindel  and  Kennel,  19713*  Tbe  estimate 
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of  n^/^  as  high  as  3 x 10  2 suggests  that  striations  may  drive  F 
region  images  at  z « 300  km  with  nx  ^ 10 4 cm  3 for  X = 1.0  km  and 
solarmax  conditions.  (The  highest  values  for  ni  can  be  found  above 
the  barium  cloud  because  of  the  large  size  of  nQ  in  the  F region  con- 
ductivity peak;  however  under  twilight  conditions  one  would  expect 
more  pronounced  E region  effects.) 

The  eigenvalue  analysis  with  partial  coupling  to  the  F region 

indicates  for  E b/E  * = 4.0  that  y is  not  significantly  decreased 
P P r 

(if  at  all)  at  short  wavelengths;  on  the  other  hand,  for  E^  /E^  =0.2, 

with  X ^ 0.2  km,  Yr  (with  T = 103  °K,  £=  5 mV/m)  appears  at  least 
20$  below  its  largest  values.  The  absence  of  a short  wavelength  cut- 
off for  Y with  E b/E  1 = 4.0  indicates  the  need  for  consideration  of 
P P 

finite  Larmor  radius  effects  and  non-linear  effects  on  the  spectral 
structuring.  Further,  for  the  same  barium  cloud  at  higher  altitudes 
(and  hence  with  lower  Pedersen  conductivities)  diffusion  damping 
effects  would  tend  to  be  more  significant. 

For  finite  temperatures  Eq.  (29)  yields  the  growth  rate: 
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The  shaping  of  the  zero  temperature  spectrum  (as  a function  of  k)  is 
generally  in  agreement  with  the  intuitive  picture  of  more  coupling  to 
the  ionosphere  and  hence  more  damping  at  larger  wavelengths  (short 
circuiting  phenomena).  Indeed,  the  term  linear  in  £/Lx  represents 
the  growth  rate  for  a barium  cloud  coupled  to  an  ionosphere  with  complex 
conductivity  E *(k)  + i L^(k)  at  the  same  potential.  The  slight 
increase  in  damping  (or  decrease  in  growth)  between  the  finite  tem- 
perature growth  rate,  y^_,  compared  with  the  zero  temperature  growth 
rate,  -vir,at  short  wavelengths  can  be  attributed  to  the  term  linear 
in  k2?  and  corresponds  to  diffusion  damping  in  the  y-direction  for  a 
barium  ion  inhomogeneity  with  stationary  electrons  and  counter- streaming 
ionospheric  ions  (see  Appendix  C for  a discussion  of  this  effect.) 
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Appendix  A - Zero-order  ionospheric  density  gradients 

In  this  Appendix,  we  estimate,  within  the  framework  of  Eqs . (6 ) 


and  (7),  modifications  on  the  background  ionospheric  density  at  various 

altitudes  consistent  with  the  zero-order  barium  density  profile 

= a k (x)  and  e = E(x).  We  will  assume  a steady  state  for  the 
p po  - - 

region.  This  probably  leads  to  an  overestimate  in  inhomogeneities 
because  from  Eq.  (A-2),  below,  the  time  necessary  to  establish  an 
inhomogeneity  in  response  to  = a^(x),  E = E(x)  varies  as  (v/u^)"1 
and  hence  can  be  very  long  for  the  F-region. 

Since  V x E = 0,  6 (Vty  ) /bx  * 0 and  (7^  ) = constant.  Hence  in 

— O y O y 

the  system  moving  with  v = -c(Vt|rQ)  x &/B  for  the  region  of  interest 
Eq.  (6)  and  (7)  become: 
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and  in  particular  at  the  barium  level  one  has: 
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Here  the  barium  cloud  is  taken  to  occupy  the  region  x > 0,  E^>0, 
6ct^/&x  >0,  and  the  superscript  h refers  to  values  at  x = O',  i.e., 
just  outside  the  projection  of  the  barium  cloud  along  magnetic  field 
lines . 

ll  b ^ 

In  principle  we  can  regard  Ex  } E^  and  a^  as  known  in(A-3  ) thereby 

determining  E and  [through  (A-3)]  (x, z).  The  case  for  E ^ = 0 is 

x p X 

particularly  simple.  On  restricting  considerations  to  the  region 
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For  L (x)  a minimum,  where  one  expects  mode  growth  to  be  localized, 
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Thus  there  are  indications  of  possible  contributions  to  growth 
from  the  zero-order  F region;  however,  as  noted,  the  time  to  set  up 
a steady  state  in  the  F region  would  appear  to  be  much  greater  than 
that  for  the  E region  and  makes  these  estimates  somewhat  more  tenuous 
for  the  F region.  For  the  E region  on  the  other  hand,  the  indication 
is  that  inhomogeneities  should  be  relatively  unimportant. 
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Appendix  B - Numerical  Method 


Here,  we  briefly  outline  the  numerical  procedure  called  BODEL 
(Dr.  E.  Graham,  private  communication)  which  was  used  to  solve  Eq.  (22). 
First,  Eq.  (22)  was  written  as  four  equations,  each  first  order  in  the 
spatial  derivative,  by  introducing  the  variables 
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so  that  we  also  have  from  Eq.  (17) 
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where  (B-l)  and  (B-2)  are  really  four  equations  since  V and  K have  real 
and  imaginary  parts  (as  does  T).  The  real  and  imaginary  part  of 
the  eigenvalue  Yi  were  taken  to  satisfy  the  two  equations 


(B-3) 


We  then  have  a sixth  order  system  with  six  equations,  six  unknowns  and 
six  boundary  conditions  (Eqs.  (23a),  (23b)  and  (23c)).  The  finite 
difference  equations  are  then  obtained  by  making  the  substitutions 


bf 


ST*  <W  £i>/4z>  £<i  <fi+i  + £i> 


at  the  i th  gridpoint  (where  f represents  one  of  the  six  variables), 


The  system  of  finite-difference  equations  was  solved  by  the  Newton- 

Raphson  iteration  method.  The  ith  of  I equations  (I  = 6),  at  the 

kth  gridpoint,  is  written  as  F (u.  ) = 0,  where  there  are  I variables 

i J ) R 

Uj,  and  M meshpoints.  Let  our  nth  approximation  to  the  solution  be 

the  I by  M matrix  u”  . . Then  for  u?*1  = u^  + 6u.  , , substitution 

j jR  J >R  J >R  j .>  R 

in  the  system  and  linearizing  gives 

Fi<uj,k>  +X)(£,Fi/e'uj,<<)6uj,k = 0 

which  is  a system  of  I equations  at  each  meshpoint.  This  I by  M 
system  is  solved  by  Gaussian  elimination  using  pivotal  condensation 
with  back  substitution.  The  matrix  is  block  bi-diagonal,  since  at 
the  k-1/2  th  point,  the  functions  depend  only  on  the  values  at  the  k 
and  k-1  8t  meshpoints.  The  Jacobian  matrix  itself,  bF/6u,  is  found 
through  numerical  differencing.  The  subprogram  which  perform  these 
operations  are  all  part  of  the  BODEL  language.  We  iterate  until  the 
maximum  residual  (value  of  F^)  is  below  an  error  criterion,  taken 
throughout  this  work  as  10_8.  The  boundary  conditions  are  treated 
similarly.  No  differencing  need  be  done  on  the  boundaries. 
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Appendix  C - The  diffusion  coefficient  for  barium  ions  in  a 
homogeneous  ionosphere  with  stationary  electrons 
We  consider  a barium  layer  with  small  inhomogeneity  in  the  x 
direction  and  an  ionospheric  layer  at  a different  altitude  but  with 
the  same  generalized  potential,  ijr(x,y).  For  the  barium  ions  we  have 


"2b 


0)  , 
— cb 


V V 

iA 


(c-i) 


where  is  the  electric  field  at  the  barium  level  and  = cu^  b/|b| . 
For  the  ionospheric  ions,  we  have 


_ v 
-i 


u>  . 
— ci 


V v 

i-i 


_!!l 

m. 


(C-2) 


where  is  the  electric  field  at  the  ionospheric  level.  Since 
(^jL^b  = ^^i  we  ^ave  on  assuming  equal  ion  and  electron  temperatures 


E = 


*Lb 


S-.i 


(c-3) 


from  which  on  integrating  over  the  barium  and  ionospheric  layers 
separately  we  obtain 


JniVz 


i 

2 E 

-JE_ 
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(C-4) 


and 
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b 


b 

r E E 

Jvtd2'-V 


(c-5) 


Since  the  electrons  are  stationary 

/niVz  - -/”*>dz 


and 


From  the  continuity  equation  we  have 


& 


= 0 


Since  the  spatial  inhomogeneity  is  small,  this  results  in: 


and  this  has  the  form 


kl 


-sr/V2  = Deff y2/  "bdz 


with  an  effective  diffusion  coefficient,  ®e££> 


eff  eB 


L1 

_E 

s’3  + 

P P 


vb/ucb 


t+t 


vb7“^r 


(c-6) 


The  result  of  Eq.  (31)  follows  from  using  = E^(k)  + iE^(k)  and 
taking  the  real  part  of  De^. 
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Appendix  D - Partial  coupling  to  the  F region 

From  Eq.  (13)  on  assuming  no  current  as  z -»  one  obtains 

* 

Z Z 

t -i-  dz'/"  r (D-d 

•00  No  •/CD 

On  writing 


r = r + r 


where 


r 

o 


and 


r 

1 


ikL. 


v/coc/ll+ 


Yi 


ikL,  v/co 

C 


we  may  expand 


♦,  ■ tio  + ♦ 


ll 


(D-2) 
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to  obtain 


z z 

♦io(z)  - W00)  -/  dz' f r WOdz' 

./»  lln  J°  ° 


♦n(«)  - txiC00)  =f  — - — dz* f T i|f11(z// )dz" 

iio 


z z 

*/  dz*  f Tx^0{z)dz 

00  ||n  J» 


where 


♦i(*x)  = ♦io(*i)  = 1 , 

♦n(zi ) = 0 

and  Z!  denotes  the  altitude  of  the  "top"  of  the  barium  cloud. 

If,  as  is  expected  for  partial  coupling  to  the  F region,  one 
♦io(®)  ^ 0,  then  ir10(z)  monotonlcally  decreases  in  magnitude  as  z 
increases.  We  may  use  this  to  show  from  (D-j5a)  that: 

+io(°°)  * 1/Cl  + Qf(zi)]  , 


where: 


kk 


(D-3a) 


(D-3b) 


has 


(D-U) 


(D-4a) 


V 


By  similar  algebra,  for  any  z*  ■Zz^  < Z*  < 00  we  may  show 

♦loO  * fl  -f  dz'[  F dz")  /A  (D-5) 

V J Z*  a||0  J 7*  ° / 

where 

dz 

The  right  hand  side  (D-5)  becomes  the  right  hand  side  of  (D-4)  for: 
rz  i zi 

ih/  -;i— d*7  r dz"  « 1 (D-6) 

A*  llo  A* 

and 

c*(z*)  « 1 . (D-7) 

On  the  other  hand  for  the  transition  from  partial  to  total  coupling 
one  seeks  by  virtue  of  (D-4)  that: 

<*(zl ) ~ 1 • (D-8) 

We  identify  z*  with  the  bottom  of  the  F region  conductivity 
maximum  and  assume  negligible  potential  drop  between  z*  and  °°.  [This 
is  equivalent  to  condition  (D-7).]  For  z * = 280  km  and  our  solarmax 
conditions  the  left  hand  sides  of  (D-6),  (D-7)  and  (D-8)  are  respectively 


z* 

'/  ‘ 


+ at  ( z* )]  . 
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0.04/X2,  O.O5/X2  and  1/X2,  X in  km.  For  z*  = 270  km  and  our  solarmin 
conditions  the  left  hand  sides  of  (D -6),  (D-7)  and  (D-8)  are  respectively 
5 x 10  3/X2,  6 x 10  3/X2  and  3.2  x 10  2/X2.  Hence  on  neglect  of  terms 
of  relative  size  ^ .2, 

♦io(®)  = 1/C  1 +Q'(z1)]  (D-9a) 


for  X ^ 0.5km  and  solarmax  conditions  or  X ^ 0.2  km  and  solarmin 
conditions.  Since  (D-9a)  results  from  replacing  ♦10  (O  by  ty10  (°° ) 
within  (D-3a)  one  expects  to  within  similar  accuracy: 

♦10(2)  = " l"T"ofz~)'  * Zl  < 2 < “ * (D-9b) 


We  next  demonstrate  that  in  these  wavelength  ranges  the  contri- 
butions from  'l'n(z)  are  small.  The  imaginary  part  of  is  bounded 
in  magnitude  by  with  6 = - kLx  v ir,  while  the  real  part 

of  r3  is  bounded  by  r^6/2.  To  simplify  the  analysis  we  consider  ra 
to  be  approximated  by  its  imaginary  part. 

On  replacing  ( z"  ) by  'I'nC00)  in  (D-3b)  and  defining 


d.y  rxt10(2#  )dz" 


(D-10) 


we  obtain 


*11 (®) 


P(*l) 

TToTzTT 


(D-ll ) 


hC 


which  itself  can  be  used  to  obtain 


♦ii(*)  - - P(*i)  + p(z) 

This  procedure  is  self-consistent  if  i]rn(z)  ^ constant  for 
z > z*  and: 


00 

/* 

r # , , 

/ rQ+ll (“)dz 

> 

J r0+n(z  )dz 
Z1 

From  (D-ll)  and  (D-12)  we  have: 


C +11  (z  )-+n  (°° )]  /+n  (°° ) = a(z)-p(z)/P(z!) 

- P(z)a(z1)/p(z1). 

Since  ^1+10(2)/!^  is  monotonically  decreasing  with  increasing 
by  successive  applications  of  the  mean  value  theorem,  it  can  be 
that  P(z)/P(zi)  < q'(z)/q'(z1).  Hence  on  using  the  estimates  for 
(d-6)  - (D-8)  for  z > z*: 

|Uii(z)  - +11  (*)] /'I'll (°° )|  < 0.2. 


(D-12) 


(D-13) 


(D-lUa) 


shown 


Further,  from,  (D-ll)  and  (D-12)  one  has: 


r 00  r Z* 

For  solarmax  conditions  we  take  or(zi)  ^4  and  / T dz'/l  T dz*  ^ 30> 

“'z*  ° z m.  ° 

oo  z* 

for  solarmin  conditions  we  take  a'(zi)  £ 1 and/  T dz' / I T dz'  f**  6: 

Jz*  ° 4 ° 

hence  (D-13)  is  satisfied. 

To  show  that  our  basic  expansion  (D-3a)  and  (D-3b)  is  valid  we 
observe  from  (D-12)  and  (D-9b): 

-i  *n(z)  _ _i0(z)  , *P(zi) 

r z ” [lWzi )]  i|fxo(z)  l-to(z)  + "TioTzTT  * D-15) 

ip(zi ) 

Further  we  have  r(zx)  = 0 , r(®)  = j — ^ < 0,  and 


T'(z)  = (-ip')  + C»(-1P)  (p'/P-a'/tt) 

(1«)2 

with  -ip  > 0 and  - iP  * <0.  Since  v^’lrio/u)c^  is  a monotonically 
decreasing  function  for  increasing  z , by  using  the  mean  value  theorem 
on  pVp  one  can  show 


p'/P  - ot'/a  < 0 . 

Hence^r^z)  <0  and  |r(“)|  > |r(z)|,  z < <*>,  or  equivalently: 
I ♦n(z)/tio(*)|<  I ^ii(c°)/,l,io(00)| 


for  z < «>.  Thus  the  expansion  (D-2)  is  valid  to  the  extent  that 

I r(°° )|  = I - P(zi)|  « 1. 
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But  we  note  on  decomposing  - P(z!)  into  component  double  integrals 


Txdz" 


(D-l 6) 


The  first  term  is  approximated  by 

Of  ( z*  )kLx  < > 

1 1 + aj  Z;l) 


(D-17) 


where  < > represents  an  average  over  the  region  z > z*. 

The  second  term  is  approximated  by 


<*(  zi) 

' 1 1 + a(z!) 


kLx  < vi/u)ci  > 


'ir 


(D-18) 


The  third  term  is  approximated  by 


- i I kLx  [vi/U)ci(200)]/Ylr  , (D-19) 

where  I is  defined  in  (D-6)  and  [ (200 )]  denotes  the  value  of 

v./w  . at  z = 200  km. 
i ci 

Expressions  (D-17)  - (D-19)  are  peaked  at  large  k.  Under  solar- 

maximum  conditions  with  X ^ 0.5  km,  we  have  < v./w  , > < 10~2. 

i ci  ’ 

kL^  ^ 12.6,  ot  (z*)  ^ 0.2,  or  (zx)  £2*.,  Cvi/cwci(200)]  = 0.03,  I £ 0.l6 
and  vlr  w - 0.9.  Then  (D-17),  (D-l8)  and  (D-19)  are  respectively 
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bounded  by  0.01  i,  0.12  i,  and  0.07  i,  and  |r(«)|<0.22.  Under 
solarminimum  conditions  with  X ^ 0.2  km,  we  have  < v. /uj  J > < 2 x 10~3, 
kLx  ^ 31.4,  a(z*)  ^0.15,  a(Zl)  < 0.81,  [vi/u)ci(200)]  = 0.014  and 
I * 0.125.  For  L b/£  1 = 4.0  (and  vlr  « - 0.9),  (D-17),  (D-l8)  and 
(D-19)  are  respectively  bounded  by  O.Oli,  0.04i  and  0.06i  and 
I r(®)|  <0.12.  For  E b/E  * = 0.2  (and  yir  «=*  - 1/3)  we  have  |r(°°)|  <0.33* 
Hence  lr(°°)|  < 1 in  all  cases  and  the  expansion  (D-2)  with  i|i10(z) 
given  by  (D-9b)  appears  valid  over  the  asserted  wavelength  ranges. 
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Appendix  E - Coupling  from  counter- streaming  ion  diffusion  for 


short  wavelengths 

For  sufficiently  short  wavelengths  the  background  ionosphere  is 
effectively  homogeneous.  Here  we  provide  an  estimate  of  the  effect 


a q are  the  local  ionospheric  Pedersen  conductivity  and  barium 
Pedersen  conductivity.  A somewhat  similar  estimate  for  small  barium 
clouds  has  been  carried  out  by  Perkins  et  al.  [19733*  However,  imaging 
effects  involving  kLxv/ouc  were  neglected  and  the  assumption  was  made 
that  the  region  of  ionospheric  homogeneity  was  sufficiently  large 
that  its  integrated  Pedersen  conductivity  could  be  comparable  with 
the  integrated  barium  Pedersen  conductivity. 

The  first  order  potential  is  described  by  Eq.  (13): 


of  coupling  for  a barium  cloud  such  that  a*  /a*5  « 1 where  a*  and 


(E-l) 


the  barium  cloud  ( z.±  > z > zq)  one  has: 


r = rt  = ri  + k2apQ(Yi  + 1)/(Y;l  - 


a 
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where  and  n^  are  the  barium  number  density  and  the  ionospheric 
number  density  respectively.  Within  each  of  the  three  regions  by 
assumption  of  homogeneity  one  has 

b2'!’! 

■ °llo  -5P~  • 


Across  Z!  and  zQ  one  has  that  fx  and  btyx/bz  are  continuous. 

Further  one  has  the  requirements  ^(z  = zx  + Ljj  )/'l'i(z1)  « 1 and 
i|ix(z  = zq  - L||  ) ( z0 ) <<:  1 where  L j is  an  ionospheric  scale  length; 
this  leads  to  the  variation 


“(ri/an  )l/2(z-zi) 

~ e ^lo  ' , z > zx 

(E-2) 

(r±/cf,,i  )l/2(z-z Q) 

i|rx  ~ e Ho  , z < zq 

(E-3) 

with  Re(T  /a  * )x/2  > 0. 

1 llo 

From  (E-l ) we  obtain 

Zi+L|| 

Ttidz  = 0 


On  anticipating  that  ii(z)  « constant  for  z 
verified  self-consistently  for  Yi^-l  and  o 
range  of  k such  that  Ref^i/c  ^ )“i/s  < 


< z < Zx  (as  may  be 

k » a*-)  for  a wide 
P P 

we  obtain,  with 


> 
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rtL|,  „ + )^lz=zi  ♦ rydyc,*  = o 

(E-4) 

It  is  clear  for  order  of  magnitude  estimates  that  we  can  evaluate 
each  of  the  ionospheric  contributions  at  the  barium  level. 

On  assuming  ap^/ap^  « 1 (and  consequently  Yi  ^ -l),  Eq.  (E-4) 
results  in: 
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We  note:  y/(c£/BLx)  = Yi  ^-1  + 2i 

one  has  2Tk/e£»  1 and 


For  short  wavelengths 


2T 


lc 


a 

PQ 

b 

a 

po 


kLx(vb/tt,cb"v/<Bc) 


Ll  + ('kLj_v/a)c)^J 


/ i \ 1/2 

Jg  (°\\o  \ 1 

kLHb  y apJ  y Ll  + (kLjV/u^W* 


^(kLxvb/(ucb)|Cl  + ( kL^v  /(D^  )2] 1 / 2 + ll3-/2 
- j[l  + (kL^v/u^)2]1/2  - l}1/2) 

As  k increases  for  either  kL  v/oi  « 1 or  kL  v/co  » 1,  y 

x c x c ' r 

decreases.  Hence  the  maximum  value  of  Yr  for  short  wavelengths 
occurs  for  k such  that  the  solution  is  just  localized.  This  is 
estimated  by  Re(0||^/r^  = L|j  with  Yi  = -i  and  yields 

k = uj^/l^v^Ti2 

Li|  / v /a>  \ 1/2 

with  T|  = J§  " [ e . Ce  ) . Characteristically  one  has 

Lx  V Vi/luci  / 

Ve/tuce  = 10"3  Vi/uci'  Lll/Lx  = 20,71  = and  vb/ucb  = ^ Vi/uci* 
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This  results  in 
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Fig.  1 —Zero  order  Pedersen  conductivity,  Opo,  vs  altitude  at  nighttime  solarmax  and  solarmin  conditions. 
These  are  basically  the  values  given  by  Hanson  [1961  ],  correcting  the  error  in  his  curve  for  solarmin  below 
the  E region  peak.  The  curves  are  consistently  based  on  Hanson’s  [1961]  curves  for  the  ambient  electron 
density  and  ion-neutral  collision  frequency  vs  altitude  using  B = 0.4  gauss  and  an  ionic  mass  equal  to  the 
neutral  atmospheric  mass  as  given  by  Johnson  [1961].  To  obtain  the  units  of  mhos/m  multiply  the  values 
given  in  this  figure  by  10^. 
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Fig.  2 — Zero  order  parallel  conductivity,  a||0,  vs  altitude,  at  nighttime  solarmax  and  solarmin  conditions. 
To  obtain  the  units  of  mhos/m  multiply  the  values  given  in  this  figure  by  10^. 
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Fig.  3 — Ratio  of  ion-neutral  collision  frequency  to  ion  cyclotron  frequency,  vs  altitude  for  solarmax  and  solarmin  conditions. 

The  values  of  are  taken  from  Hanson  [1961]  and  is  computed  using  B = 0.4  gauss  and  an  ionic  mass  equal  to  the  neutral  atmo- 
spheric mass  as  given  by  Johnson  [1961]. 


Fig.  4 — Phase  variation  in  modal  structure-zero  phase  locations  drawn  as  a function  of  altitude,  z,  given  that  at  z = 190  km 
and  y * 0 is  at  zero  phase.  Contours  are  reflected  around  y = 0 for  Bz  < 0.  a)  Case  I,  X = 0.8  km,  b)  Case  II,  X = 0.4  km, 
where  we  have  taken  £ = 5mV/m  and  T = 10^  ^K. 


CASE  I 


<c> 


Fig.  5 — Case  I,  a)  M'^(z)l  as  a function  of  altitude  and  wavelength,  X (perpendicular  to  the  magnetic  field  in  the  y direction  and 
equal  to  27r/k)  normalized  to  h^j(190)l  « 1,  b)  lviz(z)l  as  a function  of  altitude  and  wavelength,  X,  with  i//i(190)  e fi/k,  c)  (nj(z)l 
as  a function  of  altitude  and  wavelength,  X,  with  ^(190)  * &/k.  For  b)  and  c),  if  ^ “ a&/k  then  n^  and  Vjz  in  the  figure  should 
be  multiplied  by  Or. 


Fig.  7 — Case  III,  a)  \\p^(z)  \ as  a function  of  altitude  and  wavelength,  X (perpendicular  to  the  magnetic  field  in  the  y direction  and 
equal  to  27r/k)  normalized  to  h//j(liH))l  * 1,  b)  lv^z(z)l  as  a function  of  altitude  and  wavelength,  X,  with  ^(190)  *=  S/k,  c)  In^z)! 
as  a function  of  altitude  and  wavelength,  X,  with  ^i(190)  85  £/k.  For  b)  and  c),  if  \) '/^  = a£/k  then  n^  and  vlz  in  the  figure  should 
be  multiplied  by  a. 


Fig.  8 — Case  IV,  a)  JV'i(z) I as  a function  of  altitude  and  wavelength,  X (perpendicular  to  the  magnetic  field  in  the  y direction  and 
equal  to  27r/k)  normalized  to  l^^(190)l  * l,b)  lvj2(z)l  as  a function  of  altitude  and  wavelength,  X,  with  \I/U190)  *=  &/k,  c)  tai(z)l 
as  a function  of  altitude  and  wavelength,  X,  with  i^i(190)  = &/k.  For  b)  and  c),  if  « ttfi/k  then  n^ana  in  the  figure  should 
be  multiplied  by  a. 
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Fig.  9 — Case  V,  a)  hp  j(z)l  as  a function  of  altitude  and  wavelength,  X (perpendicular  to  the  magnetic  field  in  the  y direction  and 
equal  to  2n/k)  normalized  to  10^(190) I = 1,  b)  lv-iz(z)l  as  a function  of  altitude  and  wavelength,  X,  with  ^^(190)  * fi/k,  c)  kij(z)l 
as  a function  of  altitude  and  wavelength,  X,  with  \p^(  190)  * £/k.  For  b)  and  c),  if  ip^  = aS/k  then  n^  and  v^z  in  the  Figure  should 
be  multiplied  by  a. 


05 

00 


(a) 


0.5  t 


il 


^ .4 X » 02  Km 

3 ^6 


-025  J 


0.5  t 

1 


r0.00 


\ 


9 1 


-3  08  Km 


>9 — I 


0~1 — -2  2 0 Km 

4 0 Km  Q — 2^. 

0.5  12 


(b) 


-M-^4-1 1 1 

X 

\ X • 0 2 Km 

! X 


-3- 


<!f— 3-4 1 1 


\ 


4 04  Km 


-21 


-0— 2<-^~ 


\08 


-21  9 

T V^lcm/secJ/IO5 

1 2 3 

0^~ 1 1 1 

2\ 

3 2.0  Km 

+ \» 

^9 

-2^ 


Fig.  10  — Hodographs  for  Case  I.  Digits  m along  curves  indicate  values  at  altitudes  given  by  (100  + 10m)  km. 
a)  ^ii(z)  vs  ^ir(z)*  is  given  by  distance  from  origin.  The  striation  structure  has  period  X along  y. 

Locations  in  phase  with  0^(190)  at  y = 0 are  given  by  y = -£(z)X/27T  where  £(z)  is  the  polar  angle,  b)  vlzj(z) 
vs  vij^fz).  fv^2(z) I is  given  by  distance  from  origin.  Locations  in  phase  with  ^(190)  are  given  by  y = -£(z) 
X/2tt  where  £(z)  is  the  polar  angle,  c)  njj(z)  vs  n^r(z).  In^(z)l  is  given  by  distance  from  origin.  Locations  in 
phase  with  0^(190)  are  given  by  y = -£(z)X/27T,  where  £(z)  is  the  polar  angle. 
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Fig.  11  — Hodographs  for  Case  II.  Digits  m along  curves  indicate  values  at  altitudes  given  by  (100  + 10m)  km.  a)  ^]j(z)  vs  ^jr(z). 
Ii//^(z)  I is  given  by  distance  from  origin.  The  striation  structure  has  period  X along  y.  Locations  in  phase  with  ^(lOO)  at  y * 0 are 
given  by  y « -£(z)X/27 r where  £(z)  is  the  polar  angle,  b)  vlzi(z)  vs  vlzr(z).  Ivlz(z)l  is  given  by  distance  from  origin.  Locations  in 
phase  with  ^(lOO)  are  given  by  y = -£(z)\/27T  where  £(z)  is  the  polar  angle,  c)  n|j(z)  vs  n^r(z).  (n^(z) I is  given  by  distance  from 
origin.  Locations  in  phase  with  \|/^(190)  are  given  by  y * -£(z)X/27T,  where  £(z)  is  the  polar  angle. 
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Fig.  12  — Hodographs  for  Case  III.  Digits  m along  curves  indicate  values  at  altitudes  given  by  (100  + 
10m)  km.  njj(z)  vs  nlr(z).  Inj(z)l  is  given  by  distance  from  origin.  Locations  in  phase  with  1^(190) 
are  given  by  y = -£(z)X/2tt,  where  £(z)  is  the  polar  angle,  but  note  d£/dz  = 0. 
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Fig.  13  — Hodographs  for  Case  IV.  Digits  m along  curves  indicate  values  at  altitude  given  by  (100  + 10m)km. 
a)  ^if(z)  vs  ^ ir(z).  h//j(z)l  is  given  by  distance  from  origin.  The  striation  structure  has  period  X along  y. 
Locations  in  phase  with  \t^^(190)  at  y * 0 are  given  by  y * -£(z)X/27T  where  £(z)  is  the  polar  angle,  b)  vi2j(z) 
vs  v^zr(z).  klz(z)l  is  given  by  distance  from  origin.  Locations  in  phase  with  i//^(190)  are  given  by  y = -|(z) 
X/27T  where  £(z)  is  the  polar  angle,  c)  n-jj(z)  vs  nir(z).  bj(z)l  is  given  by  distance  from  origin.  Locations  in 
phase  with  i^j(190)  are  given  by  y * -£(z)X/27T,  where  £(z)  is  the  polar  angle. 
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Fig.  14  — Hodographs  for  Case  V.  Digits  m along  curves  indicate  values  at  altitudes  given  by  (100  + 10m)  km.  a)  0jj(z) 
vs  ^r(z).  h//j(z)l  i*  given  by  distance  from  origin.  The  striation  structure  has  period  X along  y.  Locations  in  phase  with 
^(190)  at  y = 0 are  given  by  y ■ -£(z)X/2tt  where  £(z)  is  the  polar  angle,  b)  vlzi(z)  vs  vlzr(z).  lvlz(z)l  is  given  by  dis- 
tance from  origin.  Locations  in  phase  with  ^^(190)  are  given  by  y * -£(z)X/ 2n  where  £(z)  is  the  polar  angle,  c)  nij(z) 
vs  nlr(z).  Ini(z)l  ia  given  by  distance  from  origin.  Locations  in  phase  with  1^(190)  are  given  by  y ■ -£(z)X/27T,  where 
£(z)  is  the  polar  angle. 
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Fig.  15  — Eigenvalues  for  Cases  (I)  - (V).  a)  the  zero  temperature  growth  rate,  -7 lr,  *n  units  of  c£/BLj_. 
"Ill”  refers  to  = 0.2,  but  with  neglect  of  image  density  transport;  b)  -7jj,  the  zero  temperature 

frequency  in  units  of  c&/BLjj  c)  7r,  the  growth  rate  appropriate  to  T * 10^  °K,  & ■ 5mV/m,  in  units  of 
cfi/BLj^.  [Eq.  (15)  with  d7^/dz  = 0,  relates  7 to  7 j - ] 


73 


74 


DISTRIBUTION  LIST 


DIRECTOR 

Defense  Advanced  Rsch  Proj  Agency 
Architect  Building 
1400  Wilson  Blvd. 

Arlington,  VA  22209 

ATTN:  Strategic  Tech  Office 

Defense  Communication  Engineer  Center 
i860  Wiehle  Avenue 
Res ton,  VA  22090 

ATTN:  CODE  R820  R.  L.  Crawford 

Defense  Communication  Engineer  Center 
i960  Wiehle  Avenue 
Reston,  VA  22090 

ATTN:  CODE  R4l0  W.  D.  Dehart 

DIRECTOR 

Defense  Communications  Agency 
Washington,  D.  C.  203  05 
ATTN:  CODE  960 

DIRECTOR 

Defense  Communications  Agency 
Washington,  D.  C.  20305 
ATTN:  CODE  480 


12  copies  (if  open  publication) 
2 copies  (if  otherwise) 

DIRECTOR 

Defense  Intelligence  Agency 
Washington,  D.  C.  20301 
ATTN:  W.  Wittig  DC-7D 

DIRECTOR 

Defense  Intelligence  Agency 
Washington,  D.  C.  20301 
ATTN:  DT-1B 

DIRECTOR 

Defense  Nuclear  Agency 
Washington,  D.  C.  20305 
ATTN:  STSI  Archives 

DIRECTOR 

Defense  Nuclear  Agency 
Washington,  D.  C.  20305 
ATTN:  STVL 


Defense  Documentation  Center 
Cameron  Station 
Alexandria , VA  223 
ATTN:  TC 


1 


T 


DISTRIBUTION  LIST 
(Continued) 


DIRECTOR 

Defense  Nuclear  Agency 

Washington,  D.  C.  20305  2 copies 

ATTN:  STT'L  Tech  Library 


DIRECTOR 

Defense  Nuclear  Agency 
Washington,  D.  C.  20305 
ATTN:  DDST 

DIRECTOR 

Defense  Nuclear  Agency 
Washington,  D.  C.  2Q305 
ATTN:  RAAE 

DIR  OF  DEFENSE  RSCH  & ENGINEERING 
Washington,  D.  C.  20301 

ATTN:  DD/S&SS  John  B.  Walsh 

DIR  OF  DEFENSE  RSCH  & ENGINEERING 
Washington,  D.  C.  20301 

ATTN:  OAD/EPS  LTC  W.  A.  Whitaker 

COMMANDER 
Field  Command 
Defense  Nuclear  Agency 
Kirtland  APB,  NM  87115 
ATTN:  FCFR 


COMMANDER 
Field  Command 
Defense  Nuclear  Agency 
Kirtland  AFB,  NM  87II5 

ATTN:  FCFR  COL.  John  P.  Hill 

Interservice  Nuclear  Weapons  School 
Kirtland  AFB,  NM  87II5 
ATTN:  Document  Control 


DIRECTOR 

Joint  Strat  TGT  Planning  Staff  Jcs 
Offutt  AFB 
Omaha , NM  68113 
ATTN:  JLTW-2 


DIRECTOR 

Joint  Strat  TGT  Planning  Staff  Jcs 
Offutt  AFB 
Omaha,  NB  68113 

ATTIM:  JPST  G.  D.  Burton 


2 


DISTRIBUTION  LIST 
(Continued) 


DIRECTOR 

Joint  Strat  TGT  Planning  Staff  Fes 
Offutt  AFB 
Omaha,  NB  68113 
ATTN:  JPST  MAJ.  J.  S.  Green 

CHIEF 

Livermore  Division  FLD  Command  DNA 
Lawrence  Livermore  Laboratory 
P.  0.  Box  808 
Livermore , CA  9^550 
ATTN : FCHUL 

COMMANDER 

National  Military  Comd  Sys  Support  Ctr 
Pentagon 

Washington,  D.  C.  20301 
ATTN:  B211 

COMMANDER 

National  Military  Comd  Sys  Support  Ctr 
Pentagon 

Washington,  D.  C.  20301 
ATTN:  DP  Director  for  CSFO 

DIRECTOR 

National  Security  Agency 
National  Security  Agency 
Ft.  George  G.  Meade,  Md.  20755 
ATTN:  W14  Pat  Clark 

DIRECTOR 

National  Security  Agency 
Ft.  George  G.  Meade,  Md.  20755 
ATTN : Frank  Leonard 

ojcs/j-3 

The  Pentagon 

Washington,  D.  C.  20301 

ATTN:  J-3  OPS  ANAL  BR.  COL.  Longsberry 

ojcs/j-6 

The  Pentagon 

Washington,  D.  C.  20301 
ATTN:  J-6 

DIRECTOR 

Telecommunications  & Comd  & Con  Sys 
Washington,  D.  C.  20301 

ATTN:  ASST  DIR  Info  & Space  Sys 


5 


DISTRIBUTION  LIST 
(Continued) 


DIRECTOR 

Telecommunications  & Comd  & Con  Sys 
Washington,  D.  C.  2Q301 
ATTN:  DEP  ASST.  SEC  Sys 

Weapons  Systems  Evaluation  Group 
400  Army  Navy  Drive 
Arlington,  VA  22202 
ATTN:  DOCUMENT  CONTROL 

DIRECTOR 

BMD  Advanced  Tech  Ctr. 

Huntsville  Office 
P.  0.  Box  1500 
Huntsville,  AL  35807 

ATTN:  ATC-T  Melvin  T.  Capps 

MANAGER 

BMD  Program  Office 
1300  Wilson  Blvd. 

Arlington,  VA  22209 
ATTN:  Plans  Division 

MANAGER 

BMD  Program  Office 
1300  Wilson  Blvd. 

Arlington,  VA  22209 

ATTN:  DACS-BMS  Julian  Davidson 

COMMANDER 

Harry  Diamond  Laboratories 
2800  Powder  Mill  Road 
Adelphi , MD.  20783 
ATTN:  AMXDO-NP 

COMMANDER 

TRASANA 

White  Sands  Missle  Range,  NM  88002 
ATTN:  EAB 

DIRECTOR 

U.  S.  Army  Baliistic  Research  Labs 
Aberdeen  Proving  Ground,  MD.  21003 
ATTN:  AM-CA  Franklin  E.  Niles 

U.  S.  Army  Communications  CMD 
C-E  Services  Division 
Pentagon  Rm.  2D513 
Washington,  D.  C.  20310 
ATTN : CEAD 


k 


DISTRIBUTION  LIST 
(Continued) 


COMMANDER 

U.  S.  Army  Electronics  Command 
Fort  Monmouth,  J.  J.  07705 
ATTN:  AMSEL-TL-ENV  HANS  A.  BOMKE 

COMMANDER 

U.  S.  Army  Material  Command 
5001  Eisenhower  Avenue 
Alexandria,  VA  22335 

ATTN:  AMCRD-WN-RE  JOHN  F.  CORRIGAN 

COMMANDER 

U.  S.  Army  Material  Command 
Foreign  and  Scientific  Tech  Center 
220  7th  St.  N.  E. 

Charlottesville , VA  22901 
ATTN:  P.  A.  Crowley 

COMMANDER 

U.  S.  Army  Material  Command 
Foreign  and  Scientific  Tech  Center 
220  7th  St.  N.  E. 

Charlottesville,  VA  22901 
ATTN:  R.  Jones 

COMMANDER 

U.  S.  Army  Missile  Command 
Redstone  Arsenal 
Huntsville,  AL  35809 

ATTN:  AMSMI-YTT  W.  G.  Freussel,  Jr. 

COMMANDER 

U.  S.  Army  Nuclear  Agency 
Fort  Bliss,  TX  79916 
ATTN:  CDINS-E 

COMMANDER 

U.  S.  Army  Nuclear  Agency 
Fort  Bliss,  TX  79916 
ATTN:  USANUA-W.  J.  Berbert 

CHIEF  OF  Naval  Research 
Department  of  the  Navy  I 

Arlington,  VA  22217 

ATTN:  CODE  464  Jacob  L.  Warner 

Chief  of  Naval  Research 
Department  of  the  Navy 
Arlington,  VA  22217 

ATTN:  CODE  464  Thomas  P.  Quinn 


5 


DISTRIBUTION  LIST 
(Continued) 


COMMANDER 

Naval  Air  Systems  Command 
Headquarters 

Washington,  D.  C.  21360 
ATTN:  AIR  5381 

COMMANDER 

Naval  Electronics  Systems  Command 
Naval  Electronic  Systems  CMD  HQS 
Washington,  D.  C.  203  60 

ATTN:  NAVALEX  034  T.  Barry  Hughes 

COMMANDER 

Naval  Electronics  Systems  Command 
Naval  Electronic  Systems  CMD  HQS 
Washington,  D.  C.  203 60 

ATTN:  IME  106-1  Satellite  Comm  Project  Off 

COMMANDER 

Naval  Electronics  Systems  Command 
Naval  Electronic  Systems  CMD  HQS 
Washington,  D.  C.  20360 
ATTN:  John  E.  Doncarlos 

COMMANDER 

Naval  Electronics  Systems  Command 
Naval  Electronic  Systems  CMD  HQS 
Washington,  D.  C.  203 60 
ATTN:  FME  117 

COMMANDER 

Naval  Electronics  Laboratory  Center 
San  Diego,  CA  92152 
ATTN:  William  F.  Moler 

COMMANDER 

Naval  Electronics  Laboratory  Center 
San  Diego,  CA  92152 
ATTN : CODE  2200  1 Verne  E.  Hildebrand 


COMMANDER 

Naval  Electronics  Laboratory  Center 
San  Diego,  CA  92152 
ATTN:  R.  Eastman 

COMMANDING  OFFICER 
Naval  Intelligence  Support  CTR 
4301  Suitland  Road,  Bldg.  5 
Washington,  D.  C.  20390 
ATTN:  Mr.  Dubbin  Stic  12 


6 


DISTRIBUTION  LIST 
(Continued) 


DIRECTOR 

Naval  Research  Laboratory 
Washington,  D.  C.  20375 

ATTN:  HDQ,  COMM  DIR  Bruce  Wald 

DIRECTOR 

Naval  Research  Laboratory 
Washington,  D.  C.  20375 

ATTN:  CODE  5460  Radio  Propagation  BR 

DIRECTOR 

Naval  Research  Laboratory 
Washington,  D.  C.  20375 

ATTN:  CODE  7127,  Charles  Y.  Johnson 

DIRECTOR 

Naval  Research  Laboratory 
Washington,  D.  C.  20375 

ATTN:  CODE  7701  Jack  D.  Brown 

DIRECTOR 

Naval  Research  Laboratory 

Washington,  D.  C.  20375  25  copies  (if  unclassified) 

ATTN:  CODE  7700,  Division  Superintendent  1 copy  (if  classified) 

DIRECTOR 

Naval  Research  Laboratory 

Washington,  D.  C.  20375  150  copies  (if  unclassified) 

ATTN:  CODE  7750,  Branch  head  1 copy  (if  classified) 

COMMANDER 

Naval  Space  Surveillance  System 
Dahlgren,  VA  22448 
ATTN:  CAFT.  J.  H.  Burton 

COMMANDER 

Naval  Surface  Weapons  Center 
White  Oak,  Silver  Spring  MD.  20910 
ATTN:  CODE  730  Tech.  Lib. 

COMMANDER 

Naval  Surface  Weapons  Center 
White  Oak,  Silver  Sprin,  Md.  20910 
ATTN:  CODE  1224  Navy  Nuc  Prgms  Off 

DIRECTOR 

Strategic  Systems  Project  Office 
Navy  Department 
Washington,  D.  C.  20376 
ATTN:  NSP-2141 


7 


DISTRIBUTION  LIST 
(Continued) 

COMMANDER 

adc/ad 

ENT  AFB  CO  80912 
ATl'N : ADDA 

AF  Cambridge  Rsch  Labs,  AFSC 
L.  G.  Hanscom  Field 
Bedford,  MA  01730 
ATTN:  LKB  Kenneth  S.  W.  Champion 

AF  Cambridge  Rsch  Labs,  AFSC 
L.  G.  Hanscom  Field 
Bedford,  MA  OT73O 
ATTN:  OPR  Hervey  P.  Gauvin 

AF  Cambridge  Rsch  Labs,  AFSC 
L.  G.  Hanscom  Field 
Bedford,  MA  01730 
ATTN:  OPR  James  C.  Ulwick 

AF  Weapons  Laboratory,  AFSC 
Kirtland  AFB,  NM  87II7 
ATTN:  DYT  LT  Mark  A.  Fry 

AF  Weapons  Laboratory,  AFSC 
Kirtland  AFB,  NM  87II7 
ATTN:  DIT  CAPT  Wittwer 

AF  Weapons  Laboratory,  AFSC 
Kirtland  AFB,  NM  87H7 
ATTN:  John  M.  Kamm  SAS 

AF  Weapons  Laboratory,  AFSC 
Kirtland  AFB,  NM  87II7 
ATTN:  SUL 

AF  Weapons  Laboratory,  AFSC 
Kirtland  AFB,  NM  87H7 
ATTN:  SAS  John  M.  Kamm 

AFTAC 

Patrick  AFB,  FL  32925 
ATTN:  TF  MAJ.  E.  Hines 

AFTAC 

Patrick  AFB,  FL  32925 
ATTN:  TF/CAFT.  Wiley 


8 


DISTRIBUTION  LIST 
(Continued) 


AFTAC 

Patrick  AFB,  FL  32925 
ATTN:  TN 

Air  Force  Avionics  Laboratory,  AFSC 
Wright-Ratterson  AFB,  Oil  45433 
ATTN:  AFAL  AVWE  Wade  T.  Hunt 

Assistant  Chief  of  Staff 
Studies  and  Analysis 
Headquarters,  U.  S.  Air  Force 
Washington,  D.  C.  20330 

Headquarters 

Electronic  Systems  Division,  (AFSC) 
L.  G.  Hans com  Field 
Bedford,  MA  01730 
ATTN:  XRE  LT.  Micheals 

Headquarters 

Electronic  Systems  Division,  (AFSC) 
L.  G.  Hans com  Field 
Bedford,  MA  01730 
ATTN:  LTC  J.  Morin  CDEF  XRC 

Headquarters 

Electronic  Systems  Division,  (AFSC) 
L.  G.  Hans com  Field 
Bedford,  MA  OI73O 
ATTN:  YSEV 

COMMANDER 

Foreign  Technology  Division,  AFSC 
Wright-Patterson  AFB,  OH  45433 
ATTN:  TD-BTA,  LIBRARY 

HQ  USAF/RD 

Washington,  D.  C.  20330 
ATTN:  RDQ 

COMMANDER 

Rome  Air  Development  Center,  AFSC 
Griffiss  AFB,  N.  Y.  13440 
ATTN:  EMTLD  Doc  Library 

COMMANDER  IN  CHIEF 
Strategic  Air  Command 
Offutt  AFB,  NB  68113 

ATTN:  XPFS  MAJ.  Brian  G.  Stephan 


DISTRIBUTION  LIST 
(continued) 

5441ES 

Offutt  AFB,  NB  68113 

ATTN:  RDPO  LT.  Alan  B.  Merrill 

Los  Alamos  Scientific  Laboratory 

P.  0.  Box  1663 

Los  Alamos,  NM  87544 

ATTN:  DOC  CON  for  Eric  Lindman 

Los  Alamos  Scientific  Laboratory 

P.  0.  Box  1663 

Los  Alamos,  NM  87544 

ATTN:  DOC  CON  for  H.  M.  Peek 

Los  Alamos  Scientific  Laboratory 

P.  0.  Box  1663 

Los  Alamos,  NM  87544 

ATTN:  DOC  CON  for  R.  F.  Taschek 

Sandia  Laboratories 
P.  0.  Box  58OO 
Albuquerque,  NM  87 115 

ATTN:  DOC  CON  for  A.  Dean  Thornbrough 

Sandia  Laboratories 
P.  0.  Box  5800 
Albuquerque,  NM  87H5 

ATTN:  DOC  CON  for  W.  D.  Brown 

Sandia  Laboratories 
P.  0.  Box  5800 
Albuquerque,  NM  87 115 

ATTN:  DOC  CON  for  D.  A.  Dahlgren,  ORG  1722 

Sandia  Laboratories 
P.  0.  Box  58OO 
Albuquerque,  NM  87 115 

ATTN:  DOC  CON  for  J.  P.  Martin,  ORG  1732 

University  of  California 
Lawrence  Livermore  Laboratory 
P.  0.  Box  808 
Livermore,  CA  9^550 
ATTN:  Tech  Info  Dept  L-3 

Department  of  Commerce 
National  Oceanic  & Atmospheric  Admin. 
Environmental  Research  Laboratories 
Boulder,  CO  80302 
ATTN:  Joseph  H.  Pope 


10 


DISTRIBUTION  LIST 
(continued) 


Department  of  Commerce 
National  Oceanic  & Atmospheric  Admin 
Environmental  Research  Laboratories 
Boulder , CO  80J02 
ATTN:  C.  L.  Rufenach 

Department  of  commerce 
Office  of  Telecommunications 
Instatute  for  Telecom  Science 
Boulder,  CO  80302 
ATTN:  Glenn  Falcon 

Department  of  Commerce 
Office  of  Telecommunications 
Tnstatute  for  Telecom  science 
Boulder,  CO  80302 
ATTN:  William  F.  Utlaut 

Department  of  Commerce 
Office  of  Telecommunications 
Institute  for  Telecom  Science 
Boulder,  CO  80302 
ATTN:  G.  Reed 

Department  of  Commerce 
Office  of  Telecommunications 
Institute  for  Telecom  Science 
Boulder,  CO  80302 
ATTN:  L.  A.  Berry- 

Department  of  Transportation 
Transportation  Rsch.  System  Center 
Kendall  Square 
Cambridge,  MA  02142 
ATTN:  TER  G.  Harowles 

NASA 

Goddard  Space  Flight  Center 
Greenbelt,  MD  20771 
ATTN:  CODE  750  T.  Golden 

NASA 

600  Independence  Ave. , S.  W. 
Washington,  D.  C.  20546 
ATTN:  M.  Dubin 

Aerospace  Corporation 
P.  0.  Box  92957 
Los  Angeles,  CA  90009 
ATTN:  T.  M.  Salmi 


13 


DISTRIBUTION  LIST 
(Continued) 


Aerospace  Corporation 
P.  0.  Box  92957 
Los  Angeles,  CA  90009 
ATTN:  S.  P.  Bower 

Aerospace  Corporation 
P.  0.  Box  92957 
Los  Angeles , CA  90009 
ATTN:  V.  Josephs on 

Aerospace  Corporation 
P.  0.  Box  92957 
Los  Angeles,  CA  90009 
ATTN:  SMFA  for  FWW 

Analytical  Systems  Corporation 
25  Ray  Avenue 
Burlington,  MA  01803 
ATTN:  Radio  Sciences 

Avco-Everett  Research  Laboratory,  Inc. 
2385  Reverse  Beach  Parkway 
Everett,  MA  02149 

ATTN:  Richard  M.  Patrick 

Bell  Telephone  Laboratories,  Inc. 
Mountain  Avenue 
Murray  Hill,  N.J.  07974 
ATTN:  Norman  J.  Zabusky 

Boeing  Company,  The 
P.  0.  Box  3707 
Seattle , WA  98124 
ATTN:  D.  Murray 

Boeing  Company,  The 
P.  0.  Box  3707 
Seattle,  WA  98124 
ATTN : Glen  Keister 

Brown  Engineering  Compnay,  Inc. 
Cummings  Research  Park 
Huntsville,  AL  35807 
ATTN:  David  Lambert  MS  18 

California  at  San  Diego,  Univ  of 
Building  500  Mathews  Campus 
3172  Miramar  Road 
La  Jolla,  CA  92037 
ATTN:  Henry  G.  Booker 


_ 


IP 


DISTRIBUTION  LIST 
(Continued) 


Cals pan  Corporation 
P.  0.  Box  235 
Buffalo,  N.  Y.  14221 
ATTN:  Romeo  A.  Deliberis 

Computer  Sciences  Corporation 
P.  0.  Box  530 
6565  Arlington,  VA  22046 
ATTN:  H.  Blank 

Computer  Sciences  Corporation 

P.  0.  Box  530 

6565  Arlington  Blvd. 

Falls  Church,  VA.  22046 
ATTN:  Barbara  F.  Adams 

Comsat  Laboratories 
P.  0.  Box  115 
Clarksburg,  MD.  20734 
ATTN:  R.  R.  Taur 

Cornel  University 

Department  of  Electrical  Engineering 
Ithaca,  N.  Y.  14850 
ATTN:  D.  T.  Farley,  Jr. 

ESL,  INC. 

495  Java  Drive 
Sunnyvale,  CA  93102 
ATTN:  J.  Roberts 

ESL,  Inc. 

495  Java  Drive 
Sunnyvale,  CA  93102 
ATTN:  V.  L.  Mower 

ESL,  Inc. 

495  Java  Drive 
Sunnyvale , CA  93 102 
ATTN:  James  Marshall 

ESL,  Inc. 

495  Java  Drive 
Sunnyvale,  CA  93102 
ATTN:  R.  K.  Stevens 


DISTRIBUTION  LIST 
(Continued) 


General  Electric  Company 
Temp-Center  for  Advanced  Studies 
8l6  State  Street 
Santa  Barbara,  CA  93102 
ATTN:  Don  Chandler 

General  Electric  Company 
Tempo-Center  for  Advanced  Studies 
8l6  State  Street  (P.0.  Drawer  QQ) 
Santa  Barbara,  CA  93102 
ATTN:  DASIAC 

General  Electric  Company 
P.  0.  Box  1122 
Syracuse,  N.  Y.  13201 
ATTN:  F.  A.  Reibert 

General  Electric  Company 
Tempo-Center  for  Advanced  Studies 
8l6  State  Street  (P.  0.  Drawer  QQ) 
Santa  Barbara,  CA  93102 
ATTN : Tim  Stephens 

General  Research  Corporation 
P.  0.  Box  3587 
Santa  Barbara,  CA  93105 
ATTN:  John  Ise,  Jr. 

Geophysical  Institute 
University  of  Alaska 
Fairbanks , AK  99701 

ATTN:  Technical  Library 

Geophysical  Institute 
University  of  Alaska 
Fairbanks,  AK  99701 
ATTN : Neil  Brown 

Geophysical  Institute 
University  of  Alaska 
Fairbanks,  AK  99701 
ATTN:  T.  N.  Davis 

GTE  Sylvania , Inc . 

I89  B Street 

Needham  Heights,  MA  02194 
ATTN:  Marshal  Cross 


14 


DISTRIBUTION  LIST 
(Continued) 


HRB-Singer,  Inc. 

Science  Park,  Science  Park  Road 
P.  0.  Box  60 

State  College,  PA  16801 
ATTN:  Larry  Feathers 

Illinois,  University  of 
Department  of  Electrical  Engineering 
Urbana,  IL  6l801 
ATTN:  K.  C.  Yeh 

Institute  for  Defense  Analyses 
4 00  Army-Navy  Drive 
Arlington,  VA  22202 
ATTN : Ernest  Bauer 

Institute  for  Defense  Analyses 
400  Army-Navy  Drive 
Arlington,  VA  22202 
ATTN : Hans  Wolfhard 

Institute  for  Defense  Analyses 
400  Army-Navy  Drive 
Arlington,  VA  22202 
ATTN:  J.  M.  Aein 

Institute  for  Defense  Analyses 
400  Army-Navy  Drive 
Arlington,  VA  22202 
ATTN:  Joel  Bengston 

inti  Tel  & Telegraph  Corporation 
300  Washington  Avenue 
Nutley,  N.  J.  07 110 

ATTN:  Technical  Library 

Itt  Electro-Physics  Laboratories,  Inc. 
9140  Old  Annapolis  Road 
Columbia,  MD.  21043 
ATTN:  John  M.  Kelso 

John  Hopkins  University 
Applied  Physics  Laboratory 
8621  Georgia  Avenue 
Silver  Spring,  MD  20910 
ATTN:  Document  Librarian 


DISTRIBUTION  LIST 
(continued) 

Lockheed  Missiles  & Space  Co.,  Inc. 

P.  0.  Box  504 
Sunnyvale,  CA  94088 
ATTN:  Dept  60-02 

Lockheed  Missiles  and  Spate  Company 
3251  Hanover  Street 
Palo  Alto,  CA  94304 
ATTN:  Billy  M.  McCormac,  Dept  52-11! 

Lockheed  Missiles  and  Space  Company 
323I  Hanover  Street 
Palo  Alto,  CA  94304 
ATTN:  Martin  Walt,  Dept  52-10 

Lockheed  Missiles  and  Space  Company 
3251  Hanover  Street 
Palo  Alto,  CA  94304 

ATTN:  Richard  G.  Johnson,  Dept  52-12 

M.I.T.  Lincoln  Laboratory 
P.  0.  Box  73 
Lexington,  MA  02173 
ATTN:  Mr.  Walden,  X113 

M.I.T.  Lincoln  Laboratory 
P.  0.  Box  73 
Lexington,  MA  02173 
ATTN:  D.  Clark 

M.I.T.  Lincoln  Laboratory 
P.  0.  Box  73 
Lexington,  MA  02173 
ATTN:  James  H.  Rannell  L-246 

M.I.T.  Lincoln  Laboratory 
P.  0.  Box  73 
Lexington,  MA  02173 
ATTN:  Lib  A-082  for  David  M.  Towle 

Martin  Marietta  Corporation 
Denver  Distribution 
P.  0.  Box  179 
Denver,  CO  80201 

ATTN:  Special  Projects  Program  248 

Maxwell  Laboratories,  Inc. 

9244  Balboa  Avenue 
San  Diego,  CA  92123 
ATTN:  A.  J.  Shannon 


DISTRIBUTION  LIST 
(Continued) 

Maxwell  Laboratories,  Inc. 

9244  Balboa  Avenue 
San  Diego,  CA  92123 
ATTN:  V.  Fargo 

Maxwell  Laboratories,  Inc. 

9244  Balboa  Avenue 
San  Diego,  CA  92123 
ATTN:  A.  N.  Rostocker 

McDonnell  Douglas  Corporation 
5301  Bolsa  Avenue 
Huntington  Beach,  CA  92647 
ATTN:  J.  Moule 

McDonnell  Douglass  Corporation 
5301  Bolsa  Avenue 
Huntingson  Beach  CA  92647 
ATTN:  N.  Harris 

McDonnell  Douglas  Corporation 
5301  Bolsa  Avenue 
Huntington  Beach,  CA  92647 
ATTN:  Tech  Library  Services 

Mission  Research  Corporation 
735  State  Street 
Santa  Barbara,  CA  93101 
ATTN:  Ralph  Kilb 

Mission  Research  Corporation 
735  State  Street 
Santa  Barbara,  CA  93101 
ATTN:  R.  Hendrick 

Mission  Research  Corporation 
735  State  Street 
Santa  Barbara,  CA  93101 
A'lTN:  Conrad  L.  Longjnire 

Mission  Research  Corporation 
735  State  Street 
Santa  Barbara,  CA  93101 
ATTN:  R.  E.  Rosenthal 

Mission  Research  Corporation 
735  State  Street 
Santa  Barbara,  CA  93101 
A'lTN:  R.  Bogus  cl  1 


17 


DISTRIBUTION  LIST 
(Continued) 

Mission  Research  Corporation 
735  State  Street 
Santa  Barbara,  CA  93101 
ATTN : David  Sowle 

Mission  Research  Corporation 
735  State  Street 
Santa  Barbara,  CA  93101 
ATTN : M.  Scheibe 

Mission  Research  Corporation 
735  State  Street 
Santa  Barbara , CA  93101 
ATTN:  P.  Fischer 

Mitre  Corporation,  The 
Route  62  and  Middlesex  Turnpike 
P.  0.  Box  208 
Bedford,  MA  01730 
ATTN:  S.  A.  Morin  m/s 

Mitre  Corporation,  The 
Route  62  and  Middlesex  Turnpike 
P.  0.  Box  208 
Bedford,  MA  01730 
ATTN:  Chief  Scientist  W.  Sen 

Mitre  Corporation,  The 
Route  62  and  Middlesex  Turnpike 
P.  0.  Box  208 
Bedford,  MA  017330 
ATTN : G . Harding 

Mitre  Corporation,  The 
Westgate  Research  Park 
1820  Dolley  Madison  Blvd. 

McLean,  VA  22101 

A'JTN:  Allen  Schneider 

North  Carolina  State  Univ  at  Raleigh 
North  Carolina  State  Univ  Campus 
Raleigh,  N.  C.  27507 

ATTN:  SEC  Officer  for  Walter  A.  Flood 

Pacific-Sierra  Research  Corp. 

1456  Cloverfield  Blvd. 

Santa  Monica,  CA  90404 
ATTN:  E.  C.  Field,  Jr. 


18 


DISTRIBUTION  LIST 
(Continued) 


Philco-Ford  Corporation 
Western  Development  Laboratories  Div 
3939  Fabian  Way 
Palo  Alto,  CA  9^303 
ATTN:  J.  T.  Mattingley  MS  X22 

Photometries , Inc . 

M2  Marrett  Road 
Lexington,  MA  02173 

ATTN:  Irving  L.  Kofsky 

Physical  Dynamics,  Inc. 

P.  0.  Box  1069 
Berkeley,  CA  9^701 
ATTN:  Joseph  B.  Workman 

R & D Associates 
P.  O.  Box  3580 
Santa  Monica,  CA  90^+03 
ATTN:  Robert  E.  Lelevier 

R & 1)  Associates 
P.  0.  Box  3580 
Santa  Monica,  CA  90403 
ATTN : Forest  Gilmore 

R & d Associates 
P.  0.  Box  3580 
Santa  Monica,  CA  90^03 
ATTN:  Richard  Latter 

R & jj  Associates 

P.  0.  Box  3580 

Santa  Monica,  CA  90^03 

ATTN:  William  B.  Wright,  Jr. 

Rand  Corporation,  The 
1700  Main  Street 
Santa  Monica,  CA  90bo6 
ATTN : Cullen  Crain 

Science  Applications,  Inc. 

P.  0.  Box  2351 
La  Jolla,  CA  92038 

ATTN:  Daniel  A.  Hamlin 

Science  Applications , Inc. 

P.  0.  Box  2351 
La  Jolla,  CA  92038 
ATTN:  D.  Sachs 


19 


DISTRIBUTION  LIST 
(Continued) 


Science  Applications,  Inc. 
P.  0.  Box  2351 
La  Jolla,  CA  92038 
ATTN:  E.  A.  Staraker 

Stanford  Research  Institute 
333  Ravenswood  Avenue 
Menlo  Park,  CA  9^025 
ATTN:  L.  L.  Cobb 

Stanford  Research  Institute 
333  Ravenswood  Avenue 
Menlo  Park,  CA  9^025 

ATTN:  Walter  G.  Chestnut 

Stanford  Research  Institute 
333  Ravenswood  Avenue 
Menlo  Park,  CA  9^025 
ATTN : Donald  Neils on 

Stanford  Research  Institute 
333  Ravensweed  Avenue 
Menlo  Park,  CA  9^025 
ATTN:  David  A.  Johnson 

Stanford  Research  Institute 
333  Ravenswood  Avenue 
Menlo  Park,  CA  9^025 
ATTN:  Charles  L.  Rino 

Stanford  Research  Institute 
333  Ravenswood  Avenue 
Menlo  Park,  CA  9^025 
ATTN:  E.  J.  Fremouw 

Stanford  Research  Institute 
306  Wynn  Drive,  N.  W. 
Huntsville,  AL  35805 
ATTN:  Dale  H.  Divis 

Tri-Cora,  Inc. 

12216  Parklawn  Place 
Rockville , MD  20852 
ATTN : Darrel  Murray 


20 


DISTRIBUTION  LIST 
(Continued) 

TRW  Systems  Group 
One  Space  Park 
Redondo  Beach,  CA  90278 
ATTN:  P.  H.  Katsos 

TRW  Systems  Group 
One  Space  Park 
Redondo  Beach,  CA  90278 
ATTN:  J.  W.  Lowry 


V 


21 


